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ABSTRACT 


Multi-frequency  modulation  (MFM)  has  been  developed  at  NPS  using  both 
quadrature-phase-shift-keyed  (QPSK)  and  quadrature-amplitude-modulated 
(QAM)  signals  with  good  bit  error  performance  at  reasonable  signal-to-noise  ra¬ 
tios.  Improved  performance  can  be  achieved  by  the  introduction  of  error  control 
coding. 

This  report  documents  a  Fortran  simulation  of  the  implementation  of  error 
control  coding  into  an  MFM  communication  link  with  additive  white  Gaussian 
noise.  Four  Reed-Solomon  codes  were  incorporated,  two  for  16-QAM  and  two 
for  32-QAM  modulation  schemes.  The  error  control  codes  used  were  modified 
from  the  conventional  Reed-Solomon  codes  in  that  one  information  symbol  was 
sacrificed  to  parity  in  order  to  use  a  simplified  decoding  algorithm  which  requires 
no  iteration  and  enhances  orror  detection  capability.  Bit  error  rates  as  a  function 
of  SNR  and  £,,/ A0  were  analyzed,  and  bit  error  performance  was  weighed  against 
reduction  in  information  rate  to  determine  the  value  of  the  codes. 
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THESIS  DISCLAIMER 


The  reader  is  cautioned  that  computer  programs  developed  in  this  research 
may  not  have  been  exercised  for  all  cases  of  interest.  While  every  effort  has  been 
made,  within  the  time  available,  to  ensure  that  the  programs  are  free  of  compu¬ 
tational  and  logic  errors,  they  cannot  be  considered  validated.  Any  application 
of  these  programs  without  additional  verification  is  at  the  risk  of  the  user. 
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I.  INTRODUCTION 


Many  modern  communication  methods  which  merge  communications  and 
computer  technology  require  efficient  and  reliable  transfer  of  digital  data  from 
an  information  source  to  a  destination.  This  requirement  is  not  always  an  easy 
one  to  achieve.  Because  this  transmission  is  through  a  physical  medium  which  is 
subject  to  various  types  of  noise,  distortion  and  interference,  the  data  output  from 
the  channel  may  differ  from  its  input.  The  demand  for  reliable  communication 
necessitates  the  detection  and  possibly  the  correction  of  those  errors  resulting 
from  impaired  transmission.  These  errors  can  be  detected  and  corrected  using  the 
principles  of  error  control  coding. 

Error  control  coding  can  be  implemented  in  a  communication  system  usually 
at  the  cost  of  data  transmission  rate.  Generali}',  the  code  that  corrects  more  er¬ 
rors  will  require  more  parity  ssmbols.  reducing  the  data  rate.  On  the  other  hand, 
a  code  which  corrects  more  errors  will  be  more  reliable  in  a  high  noise  environ¬ 
ment.  The  balance  between  data  transmission  rate  and  coding  gain  poses  an  in¬ 
teresting  problem,  one  which  will  be  addressed  in  this  thesis. 

The  subject  of  this  thesis  is  the  simulation  and  analysis  of  various  error  con¬ 
trol  codes  in  a  bandpass  Multi-Frequency  Modulation  (MFM.)  communication 
link  with  Additive  White  Gaussian  Noise  (AWGN).  MFM  is  a  signal  modulation 
format  which  is  readily  adaptable  to  many  digital  communication  link  scenarios. 
Its  efficient  use  of  the  frequency  spectrum  allows  very  high  data  transmission 
rates.  But  like  all  communication  systems.  MFM  is  subject  to  errors. 

There  are  four  error  control  codes  used  in  the  simulation,  all  of  which  are 
Reed-Solomon.  Two  error  correction  detection  capability  codes  (single  error 
correcting  double  error  detecting  (SEC  DED)  and  double  error  correcting  triple 
error  detecting  (DEC  TED))  in  combination  with  three  different  modulation 
schemes  (one  16-QAM  and  two  32-QAM  constellations)  provide  performance 
data  for  the  four  codes.  The  cffccti\ cness  of  implemented  codes,  in  terms  of  bit 
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error  rates,  coding  gains  obtained  versus  the  uncoded  link,  and  code  rates  are  the 
areas  of  study  for  this  report. 

Chapter  Two  is  a  brief  o\erview  of  MFM  as  it  pertains  to  a  16-20  kHz 
bandpass  channel,  the  communication  link  under  simulation.  Chapter  Three  is 
a  discussion  of  Reed-Solomon  error  control  coding  in  general,  including  the  en¬ 
coding  and  decoding  processes  and  other  parameters  necessary  to  define  that 
class  of  codes.  Also  included  is  a  description  of  how  the  codes  implemented  differ 
from  conventional  Reed-Solomon  codes.  Chapter  Four  describes  the  link  simu¬ 
lation  program  and  supporting  programs  and  Files,  channel  coding  and  decoding, 
how  error  control  was  implemented  into  the  MFM  environment,  and  the  results 
obtained  using  the  various  codes  compared  to  the  uncoded  channels.  Chapter 
Five  contains  the  conclusions  and  recom  lendations  of  possible  ways  to  achieve 
improx ed  performance. 


II.  THEORY  OF  MULTI-FREQUENCY  MODULATION 


A.  BACKGROUND 

Most  modulation  methods  for  bandpass  channels  hold  the  digital  signal  in¬ 
formation  in  the  air  litude  and  or  phase  of  a  earner  (e.g.,  Phase  Shift  Keying 
(PSK).  Quadrature  Amplitude  Modulation  (QA.M)).  Information  may  also  be 
transmitted  by  the  presence  of  one  of  two  carrier  frequencies  (Frequency  Shift 
Keying  ( F S K ) ) .  MFM  is  a  variation  on  these  methods.  The  signals  are  directly 
encoded,  modulated,  demodulated,  and  decoded  using  digital  signal  processing 
techniques  contained  in  the  software  of  the  host  computer. 

In  MFM.  signals  are  grouped  into  "packets"  which  are  arbitrarily  located  in 
the  frequency  spectrum,  and  consists  of  one  or  more  bauds.  An  MFM  packet  is 
shown  in  Figure  1.  on  a  display  of  frequency  and  time.  This  figure  shows  the 
nature  of  MFM:  a  structured  block  of  time  and  frequency  slots. 

This  structure  makes  MFM  a  technique  that  is  ideal!1  suited  to  allow  the  host 
computer  to  absorb  the  functions  of  modulation  and  multiplexing,  demultiplexing 
and  demodulation  [Ref.  1:  p.  3].  The  packet  consists  of  L  bauds,  each  holding 
information  on  K  tones.  The  baud  length  is  AT  during  which  K  discrete  tones 
are  transmitted.  In  MFPSK.  the  phase  of  each  tone  contains  the  information, 
and  in  MFQAM.  both  the  phase  and  magnitude  of  each  tone  contain  the  infor¬ 
mation.  These  LK  subsignals  form  an  orthogonal  signal  set.  Euch  of  these  sub¬ 
signals  may  be  independently  modulated  with  both  phase  and  amplitude 
information. 

B.  GENERATION  AND  DEMODULATION 

Generation  and  demodulation  of  MFM  is  based  on  the  symmetry  properties 
of  the  Fast  Fourier  Transform  (FFT).  That  is.  the  FFT  of  a  real-time  signal 
corresponds  to  a  conjugate  symmetric  frequency  spectrum,  and  the  inverse  FFT 
of  a  conjugate  symmetric  frequency  spectrum  corresponds  to  a  real-time  signal. 
This  means  that,  if  oniy  half  of  the  frequency  spectrum  is  used  for  in-phase  and 
quadrature  information,  and  the  other  half  is  loaded  with  its  complex  conjugate 
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Figure  1.  MFM  Signal  Packet:  (after  Ref.  1:  p.3) 


image,  then  an  inverse  FFT  of  this  spectrum  will  result  in  a  real-time  signal. 
Conversely,  the  FFT  of  this  real-time  signal  will  result  in  the  same  frequency 
spectrum.  Thus,  the  MFM  transmitter  is  a  device  which  takes  an  inverse  FFT, 
and  the  MFM  receiver  is  a  device  which  takes  a  FFT  (with  appropriate  D  A  and 
A  D  conversion). 

In  a  bandpass  channel,  the  usable  portion  of  the  frequency  spectrum  for  in¬ 
formation  is  further  reduced.  Current  MFM  applications  under  study  at  Naval 
Postgraduafe  School  deal  with  a  16-20  kHz  bandpass  channel,  with  a  D  A  and 
A  D  sampling  frequency  of  61440  Hz.  The  corresponding  parameters  for  the 
MFM  packets  are  shown  in  Table  1. 

In  the  communication  link  under  study,  the  baud  type(s)  transmitted  in  any 
give-  ,/c.ket  is  known  by  the  receiver.  Because  the  tones  containing  information 
occur  in  a  known  band  of  the  frequency  spectrum,  only  those  harmonic  numbers 
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of  the  FFT  corresponding  to  inis  band  need  be  analyzed.  This  function  is  similar 
to  bandpass  filtering,  but  noise  within  the  pass  band  is  still  present,  leading  to 
potential  errors. 


Table  1.  MFM  PACKET  PARAMETERS  (BANDPASS  16-20  KHZ) 


BAUD  TYPE 

1 

-> 

3 

4 

5 

Space  Between 
Tones  Af  (Hz) 

240 

120 

60 

30 

15 

Lower  Harmonic 
Numbcr-kl 

68 

135 

269 

537 

1073 

Lower  Tone  ( Hz) 

16320 

16200 

16140 

161  10 

1 6095 

Upper  Harmonic 
Number-k2 

S3 

1 66 

332 

664 

1 32S 

Upper  Tone  (Hz) 

1 9920 

19920 

19920 

1 9920 

1 9920 

Number  of  Tones 
With  Information 

16 

32 

64 

12S 

256 

FFT  Size 
(Baud  Length) 

512 

1024 

204S 

4096 
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III.  REED-SOLOMON  CODES 

A.  INTRODUCTION 

A  binary  block  error  control  code  with  block  length  n  which  contains  2k  code 
words  is  called  a  linear  ( n,k )  code  if  and  only  if  those  2k  code  words  form  a  k- di¬ 
mensional  subspace  of  the  vector  space  of  all  the  ^-tuples  over  the  Galois  Field 
GF(2).  A  cyclic  code  is  a  linear  block  code  in  which  every  cyclic  shift  of  a  code 
word  (left  or  right)  is  also  a  code  word.  Cyclic  codes  are  useful  because  the  en¬ 
coding  and  syndrome  calculation  processes  can  be  accomplished  easily  through 
shift  registers  with  feedback  connections  and  because  their  inherent  algebraic 
structure  makes  possible  the  use  of  several  practical  methods  for  decoding 
[Ref.  2:  p.  85].  These  codes  are  normally  expressed  as  polynomials.  A  large  class 
of  powerful  cyclic  codes  are  those  known  as  BCH  codes.  Their  importance  in 
coding  is  due  to  the  fact  that  they  are  capable  of  correcting  all  random  patterns 
of  t  errors  with  a  simple  decoding  algorithm  that  is  easily  implemented  [Ref.  3: 
P.  167]. 

Rced-Solomon  codes  form  a  special  subclass  of  non-binary  BCH  codes.  A 
conventional  (n.k)  r-error  correcting  Reed-Solomon  code  which  is  comprised  of 
elements  of  the  Galois  Field  GF(?)  has  the  following  structure: 

•  code  word  block  length  n  =  q  —  1 

9 

•  number  of  parity  symbols  n  —  k  =  2t 

•  minimum  distance  dmn  =  2t  +  1 

For  this  design,  q  is  taken  to  be  a  power  of  two  (i.e.,  q  =  2m  ).  Appendix  A  dis¬ 
plays  a  binary  representation  of  the  elements  of  GF(24)  and  GF(25). 

B.  ENCODING 

If  a  is  a  primitive  element  of  GF(2m),  then  the  generator  polynomial  of  this 
(n,k)  Reed-Solomon  code  is  given  by 
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(1) 


G(x)  =  (x  4-  a)(x  4-  a2) ...  (x  4-  a"') 

=  £o  +  £ix  +  S2X'2  +  '"  +  <?2/— i"^2f  '  +  •x‘2f 

The  generator  polynomial  roots  can  also  be  shifted  to  save  circuit  components 
[Ref.  4:  p.  174].  In  this  case,  G(x)  takes  the  form 

G(x)  =  (x  4-  g^Hx  +  (/° ~hI) ...  (x  +  oc70  +2'_l),  (2) 

where  j0  is  any  integer.  The  code  generated  by  G(x)  consists  of  all  polynomials 
with  coefficients  from  GF(2m),  of  degree  n  —  1  or  less,  and  which  are  multiples 
of  G(x).  The  message  to  be  encoded  is  a  polynomial,  U{x),  of  degree 
k  -  l  =  n  -  2t  -  1  or  less.  In  non-systematic  form,  the  corresponding  codeword 
polynomial  is 

I  (x)  =  G(x)L'(x)  =  v0  4-  v,x  4-  v2x2  4-  —  4-  vn_jx”  1  (3) 

This  is  clearly  a  multiple  of  G(x)  . 

Encoding  in  systematic  form  requires  more  computation.  To  structure  the 
codewords  with  the  information  symbols  preceeded  by  the  parity  symbols,  U(x) 
must  first  be  pre-multiplicd  by  x2'.  Then  x2l£/(x)  is  divided  by  G(x) ,  which  yields 
a  quotient  (>(x)  and  a  remainder  B{x)  of  degree  2t  —  1  or  less: 

x2,G(x)  =  (7(x)G(x)  +  B(x)  (4) 

If  £(x)  is  added  to  x2'U{x),  then  the  result  is  a  multiple  of  G(x)  which  is  the  sys¬ 
tematic  codeword  corresponding  to  C'(x)  : 
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x^r  L’(x)  +  B(x)  =  (?(x)<7(.r) 


(5) 


This  polynomial  division  can  be  implemented  in  a  shift  register  feedback  circuit 
as  shown  in  Figure  2. 


Figure  2.  Systematic  Encoding  Circuit  for  an  {n,k)  Cyclic  Code:  [Ref.  2:  p.  95J 


Using  the  circuit  in  Figure  2,  an  (n,  k)  linear  cyclic  code  with  generator 

polynomial  G(x)  as  defined  in  equation  (1)  or  equation  (2)  is  encoded  as  follows: 

1.  With  the  gate  turned  on,  the  k  information  symbols  u^,  uu  ...,  «*_,  are  shifted 
into  the  communication  channel  and  at  the  same  time  shifted  into  the  the 
encoding  circuit.  This  automatically  has  the  effect  of  pre-multiplving  U{x) 
by  x2 3' . 

2.  After  k  clock  cycles,  the  entire  message  has  been  shifted  into  the  channel  > 

and  the  n  -  k  shift  registers  contain  the  elements  of  the  remainder 
polynomial. 

3.  The  gate  is  opened  to  break  the  feedback  connections  and  the  switch  is 
shifted  to  transfer  parity  symbols. 
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4.  The  parity  symbols  are  clocked  out  onto  the  channel.  These  n  -  k  parity 
symbols  together  with  the  k  information  symbols  form  the  n-length 
codeword. 

Reed-Solomon  codes  in  systematic  form  use  the  same  encoding  circuit. 

C.  DECODING 

Decoding  Reed-Solomon  codes  can  be  performed  using  a  table-lookup 

method  or  using  iterative  algorithms  (Berlekamp,  Forney,  Euclid,  etc.)  [Refs.  4, 

\ 

5].  For  binary  BCH  codes,  these  algorithms  determine  the  location  of  errors 
(using  the  error  location  polynomial).  For  non-binary  codes,  as  in  the  case  of 
Reed-Solomon  codes,  once  the  error  location  polynomial  is  found,  the  magnitude 
of  the  error(s)  must  also  be  found  in  order  to  correct  the  received  signal  block. 

Because  every  codeword  F(v)  is  divisible  by  the  generator  polynomial  G(.r) 
and  the  roots  of  <j(x)  are  a,  a*’,  ...,  a2r,  these  are  also  the  roots  of  Wx).  As  such, 
the  following  equation  holds  for  1  <  i  <  2 1: 

l  (a'j  =  Vq  -I-  \\z‘  +  v2a2'  4-  •••  +  vn_1a("”1^  =  0  (6) 


or  in  vector  form, 


9 


1 


a 

2i 

a 

(vo.v,.v2!  =0  (7) 

a("-1,i 

By  forming  an  n  by  n  -  k  =  2r  matrix  whose  columns  arc  the  vectors  of  powers 

of  a,  /  =  1.2 . 2 r  in  the  previous  equation,  any  codeword  vector  will  be 

orthogonal  to  that  matrix.  The  matrix  thus  formed  is  the  transpose  of  the  parity 
check  matrix  for  a  conventional  (n.k)  Reed-Solomon  code.  In  other  words,  the 
parity  check  matrix  for  an  (n.k)  Rced-Solomon  code  is  given  by: 

.  2  n—  1 

1  a  a  ...  a 

1  (a')  (a  T  (^)n 

1  (a3)  (a3)2  ...  (a3r! 


I  (a2')  (a2r)2  ...  (x2l)‘ 

The  syndrome  of  a  received  polynomial  is  found  from  multiplying  the  re¬ 
ceived  polynomial  (in  vector  form)  by  the  transpose  of  the  parity  check  matrix: 


1U 


(9) 


s  =  r  •  H T  =  (v  +  e)  •  H T  =  e  •  H  T, 

where  s  is  the  syndrome  vector,  r  is  the  received  vector,  v  is  the  transmitted 
codeword  vector,  e  is  the  error  vector,  and  H  is  the  parity  check  matrix.  (Recall 
that  any  codeword  vector  v  is  orthogonal  to  HT.)  Syndrome  component  5,  can 
also  be  calculated  by  evaluating  the  received  polynomial,  R{x),  at  a',  for 
i  =  1.2.. ..,2 1.  Once  the  syndrome  components  are  found,  they  are  used  in  the  it¬ 
erative  calculations  needed  to  find  the  error  location  polynomial,  &(x),  where 

(7 (V)  =  (1  ~  /Jj A-)(l  +  /kv)  -  (1  +  M 

2  v  (1°) 

—  (7q  -t  CjA  +  -b  (7AA  "b  +  <7V.Y 

Here,  the  roots  of  a{x)  (/?,  h  /id, ...  .  /J->  )  arc  the  error  location  numbers,  and  v  is 
the  total  number  of  errors  in  the  received  codeword  (for  v  <  t).  Then  defining 


.Z(A')  =  1  -r  (Sj  -f-  (7  1  )A'  +  (.So  "b  (7 ) .*> j  +  C t )V  +  ••• 

■+•  (y.  +  +  ('2sv-2  +  +  c-v,)a-v. 


(11) 


the  error  magnitude  at  location  ft,  —  a”  is  [Ref.  2:  p.  175]: 

Z(PT') 

ei, = — - • 

M  a  +m') 

1=  1  j*l 


(12) 


This  demonstrate',  the  ease  of  implementation  of  the  iterative  algorithms:  the 
disadvantage  is  that  the  iterative  solution  does  not  provide  high-speed  decoding. 
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On  the  other  hand,  a  table-lookup  method  can  provide  high-speed  decoding,  but 
requires  too  much  memory  or  overly  complicated  logic  circuitry  even  for  moder¬ 
ate  code  lengths  [Ref.  5’  p.  1].  However,  some  modifications  to  Reed-Solomon 
codes  can  allow  a  faster  decoding  procedure. 

D.  MODIFIED  REED-SOLOMON  CODES 

There  are  several  ways  to  modify  Reed-Solomon  codes  to  fit  specific  applica¬ 
tions.  The  block  length  can  be  extended  or  shortened,  and  the  number  of  parity 
or  information  symbols  can  be  varied.  The  (2m  —  1,  2m  —2 1  —  2)  codes  utilized  in 
the  simulation  carry  fewer  information  symbols  (more  parity)  than  the  conven¬ 
tional  code-*,  reducing  the  code  rate.  The  benefits  lie  in  the  enhanced  error 
correction  detection  capability  and  speed  of  decoding. 

1.  The  (2"  -  1, 2m  -  4)  Reed-Solomon  Codes 

By  lepiacing  one  of  the  information  symbols  in  the  conventional 
(2m— 1.  2"' -3)  code  with  a  parity  symbol,  the  single-error  correcting  code  can 
also  detect  all  double  errors.  This  modified  SEC  DED  code  has  the  following 
characteristics: 

•  codeword  block  length  n  —  q  —  1 

•  number  of  parity  symbols  n  —  k  =  3 

•  minimum  distance  <4m  —  4 

Generation  of  this  code  is  similar  to  conventional  codes,  except  that  now 
the  generator  polynomial  has  an  additional  factor: 

G(.y)  =  (x  +  1  )(-Y  +  a)(.r  +  a2) 

fa3  4-  y]  ]x  +  a10x2  +  er3,  for  the  (15,12)  code  (13) 

|a3  -f  ynx  +  an.x2  +a*3,  for  the  (31.28)  code 
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The  encoding  circuit  for  these  codes  is  still  that  of  Figure  2.  The  major  difference 
between  these  modified  codes  and  the  conventional  codes  lies  in  the  decoding. 

Reference  6  details  a  fast  decoding  technique  for  an  extended  SECDED 
Reed-Solomon  code  based  solely  on  the  syndrome  components  (i.e.,  iterative  al¬ 
gorithms  are  not  required).  The  same  decoding  technique  applies  to  the  modified 
codes  used  in  the  simulation,  since  the  modified  codes  are  shortened  versions  of 
the  extended  codes.  For  the  simulation,  the  blocklength  of  the  extended  codes  is 
shortened  by  three  symbols  to  better  match  the  number  of  tones  per  baud  type, 
and  to  facilitate  encoding. 

1  he  extended  codes  of  Reference  6  are  formed  by  the  addition  of  a  three 
by  three  identity  matrix  to  the  parity  check  matrix  of  the  code  generated  by  the 
generator  polynomial  of  equation  (13).  This  parity  check  matrix  is  then  short¬ 
ened  by  deleting  the  three  right-most  columns  for  use  in  the  (2m  -1.  2m  -4)  sim¬ 
ulation  codes: 


H 


i  0  0  1  1  1  ...  1 

0  10  1a  (a)2  ...  (af~4 

0  0  1  1  a“  (a.‘)“  ...  (a ~f  4 


(14) 


The  decoding  algorithm  for  these  codes  is  as  follows: 

1.  Calculate  the  syndrome  components  (in  vector  form,  s  =  (50,  5:)).  If  all  are 

zero,  decide  no  errors  occurred. 

2.  If  a  single  error  occurs  in  one  of  the  first  three  positions  of  the  received  vec¬ 
tor  (which  corresponds  to  the  identity  matrix  positions  of  H).  then  the  syn¬ 
drome  has  only  one  non-zero  element.  The  error  has  magnitude  equal  to  the 
non-zero  syndrome  component  and  location  corresponding  to  the  syndrome 
component  number. 

3.  A  single  error  in  any  other  position  satisfies 
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*1  *2  3 

with  the  error  in  position  i  of  magnitude  50. 

4.  If  a  single  syndrome  component  is  zero,  or  if 


(15) 


(16) 


then  decide  that  at  least  two  symbol  errors  occurred. 

5.  Once  errors  are  corrected,  recalculate  the  syndrome  to  ensure  all  components 
are  zero.  If  not,  decide  that  two  or  more  symbol  errors  occurred  which 
yielded  the  same  syndrome  as  a  single  error.  (Note:  this  step  was  not  in¬ 
cluded  in  Ref.  6  but  was  found  necessary  to  detect  all  double  errors.) 


As  an  example,  suppose  that  for  the  (15.12)  code  a  single  error  of  magnitude  a9 
occurred  in  the  fifth  position  of  the  received  block  r.  Then  the  error  vector  is  given 

by  e  =  (O.O.O.O.  a-.0 . 0).  and  the  syndrome  is  s  =  r  •  Hr  =  e  •  HT  =  (a9,  a11,  a13). 

Then  y  s0  =  y  y  —  a:  =  a5-3,  corresponding  to  an  error  of  magnitude  y,  —  a0  in  the 
fifth  position. 

2.  The  '2'"  -  1 , 2"'  -  6)  Reed-Solomon  Codes 

The  (2W  -  1, 2m  -  6)  codes  arc  obtained  by  making  the  same  alteration  to 
a  conventional  double-error  correcting  code  as  was  the  case  for  the  SEC  DED 
code:  sacrifice  one  information  symbol  to  parity.  This  now  allows  the  detection 
of  all  triple  errors,  and  as  will  be  shown  later,  the  correction  of  some  triple  burst 

errors.  The  characteristics  of  the  DEC  TED  codes  are: 

•  code  word  block  length  n  =  q  —  1 

•  number  of  parity  symbols  n  —  k  =  5 

•  minimum  distance  dmin  =  6 


The  generator  polynomial  for  this  code  also  has  an  additional  factor  over 
the  conventional  code,  and  using  shifted  roots  is  given  by: 
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for  the  ( 1 5,10)  code  (17) 
for  the  (31 ,26)  code 


6(.v)  =  (x  +  a  2)(x  +  a  ')(x  4-  l)(x  4-  a)(x  4-  a2) 
f  1  4-  a4x  4  a!1x2  4  a1  'x3  4-  a4x4  4-  x5, 

=  ll+a13X  +  a17x2  +  a17x3  +  al3X4  +  X5, 


Once  again,  these  codes  are  encoded  in  systematic  form  using  the  circuit  of  Figure 
2.  Decoding  is,  however,  somewhat  more  involved  than  for  the  SEC  DED  codes. 

Reference  5  gives  a  detailed  description  and  proof  of  a  fast-decoding 
technique  for  these  codes  which  is  again  solely  based  on  the  syndrome  of  the  re¬ 
ceived  signal.  The  parity  check  matrix  of  these  codes  has  an  additional  row  over 
the  conventional  codes  and  is  given  by: 


—2  ,  -2,2 
:  (a  )  • 

"ST 

i 

1 

!  (a  ) 

■  (a  ) 

1  1 

1 

a  (x)2 

-  (»>"-' 

v.  (a  ) 

(18) 


The  steps  of  the  decoding  algorithm  are  as  follows: 

1.  Compute  the  syndrome  (in  sector  form,  s  =  (s_2,  s_,,  $0, 5,,  s2)).  If  all  compo¬ 
nents  are  zero,  decide  that  no  errors  occurred. 

2.  If  more  than  two  components  of  the  syndrome  are  zero,  decide  that  at  least 
three  symbols  are  in  error. 

3.  Compute  >•,,  y2,  and  y3,  where 


>'  1  ~  sls~2  +  s-ls0 

(19) 

>2  =  s2s-2  +  s0 

(20) 

3  =  s0sl  +  *2*-l- 

(21) 
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If  Vi  =  V:  =  >'3  =  0  then  calculate  or'  —  —  .  Correct  a  single  error  at  location 
/  of  magnitude  5C. 

4.  If  •/,.  y;.  and  y,  arc  not  all  zero  but  at  least  one  of  them  equals  zero,  decide 
that  at  least  three  symbols  are  in  error. 

5.  If  none  of  the  y/s  are  zero,  compute  k  =  c  b2  and  the  trace  of  k.  where 
b  —  y:  y,.  c  =  y3  y,,  and  the  trace  of  k.  an  element  of  GF(2m).  is  defined  as 


6. 


Tr(A-)  = 


w-l 


1=0 


(22) 


If  Tr(A')  =  1  ,  decide  that  at  least  three  symbol  errors  occurred. 

1  f  T r ( A; )  —  0 .  assume  that  a  double  symbol  error  occurred  and  solve 
v-  -f  bv  -r-r  =  0  to  find  the  roots  a’  and  a;.  Then  the  first  error  is  in  the 


/"  position  with  magnitude 


e.  = 


( \nr  4-  5. ) 


(23) 


and  the  second  error  is  in  the  j'h  position  with  magnitude 

=  T)  +  e\  ■ 


(24) 


7.  Once  errors  arc  corrected,  recalculate  the  syndrome  to  ensure  all  components 
are  zero.  If  not.  decide  that  at  least  three  symbol  errors  occurred  which 
yielded  the  same  syndrome  as  a  double  error.  (Note:  this  step  was  not  in¬ 
cluded  in  Ref.  5  but  was  found  necessary  to  detect  all  triple  errors.) 


Reference  5  does  not  examine  the  code  s  capability  of  triple  error  cor¬ 
rection.  but  the  correction  of  certain  triple  burst-errors  is  possible.  Suppose  that 
three  errors  of  magnitude  c,  occur  in  consecutive  positions  of  the  received  signal. 

the  rirst  at  position  /,  /  =  0.1 . n  —  1.  (Note  that  wrap-around  errors  can  also 

be  corrected,  e.g..  errors  in  the  (/?  —  1  ),f.  O’*  and  1!'  positions.  )  The  emor 
polynomial  £(.v)  is  of  the  form 


£(.y)  =  C|.y'  4-  Cj.v,+  1  4-  e].v'+*. 


(25) 


or  in  vector  form,  e  =  (0.0,...,0,£’,.£’1,  e,.0,...,0).  Multiplying  the  error  vector  by  the 
transpose  of  .he  parity  check  matrix  gives  a  syndrome  vector  of  the  form 


X_2 

cI(«-*W2<*,)  +  cr2,/+2>) 

s-\ 

{^(oc  1  +  +  a  (,+‘l) 

*0 

— 

ei 

*1 

^(a*  +  a!+1  +  a!+2) 

s2 

el(«2,'  +  a2^i,  +  «2,<'+2)) 

(26) 


From  these  relations,  it  is  clear  that  if  this  type  of  triple  error  occurs,  then  the 
error  locations  and  magnitude  can  be  calculated  solely  from  the  syndrome  if  the 
following  simplified  conditions  exist: 


s-2 

-  — 

CyU.  ( I  +  a  +  0(  ) 

s-l 

\l  +  a  1  +  a  2) 

s0 

= 

s  1 

CjOC  ( 1  Cl  +  OF) 

*2 

L  j 

e]a2,(l  +  a2  4-  a4) 

(27) 


Note  that  the  addition  of  elements  of  GF(2'")  vary  with  the  value  of  m,  so  the 
exact  relations  between  syndrome  components  and  error  locations  vary  with  the 
code  used.  As  an  example,  for  the  (15,10)  code  with  elements  from  GF(24),  the 
relations  will  be 
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^~2'U) 


eA*'°) 

evy2,{^) 


Then  to  check  for  this  type  of  error  in  the  (15,10)  code,  the  syndrome  components 
must  satisfy  =  (y,)-V  and  s_2s2—  (y,):a6.  The  error  magnitude  is  then  the 
value  of  Sq  and  the  starting  location  can  be  obtained  from  any  of  the  other  syn¬ 
drome  components,  say  by  dividing  st  by  %a]0  and  taking  the  exponent  of  u  . 
This  capability  allows  correction  of  an  additional  240  errors  [15  triplets  x  16  ele¬ 
ments  in  GF(24)]  for  the  (15.10)  code  and  992  errors  [31  triplets  x  32  elements  in 
GF(2')]  for  the  (31.26)  code. 


IV.  SIMULATION  AND  RESULTS 


A.  BACKGROUND 

The  MFM  link  is  simulated  in  the  Fortran-77  program  "MFMLINK".  For  a 
given  baudtype  and  code  (uncoded  16-QAM,  uncoded  32-QAM,  and  (15,10), 
(15,12),  (31,26)  and  (31.28)  Reed-Solomon  codes),  the  program  chooses  random 
data  from  GF(24)  for  16-QAM  or  from  GF(25)  for  32-QAM.  Uncoded  QPSK  was 
also  included  for  comparison  with  the  QAM  results. 

For  the  cases  with  error  control,  this  data  is  encoded  into  Reed-Solomon  code 
word  blocks.  For  each  baud,  the  elements  of  GF(2m)  are  then  transformed  into 
the  appropriate  QAM  symbol  with  each  tone  representing  one  element  of 
GF(2m),  loaded  into  the  frequency  spectrum  (with  conjugate  symmetry)  and 
transmitted  with  the  invcrse-FFT.  At  this  point  (in  the  time  domain),  white 
Gaussian  noise  is  added.  The  FFT  is  then  taken,  and  each  FFT  sample  within 
the  information  band  (from  tones  k]  to  k2)  is  transformed  into  an  element  of 
GF(2'7). 

For  the  coded  cases,  these  k2-kx+ 1  symbols  are  divided  into  code  word 
blocks,  and  each  block  is  decoded  separately.  The  output  of  the  link  is  then 
compared  to  the  input  for  error  calculation.  A  block  diagram  of  the  link  is  shown 
in  Figure  3.  A  summary  of  the  Reed-Solomon  codes  used  is  shown  in  Table  2, 
with  code  rate  R  =  kjn  and  spectral  efficiency  =  Rt]unc  (i]unc  is  the  uncoded 
MF-QAM  spectral  efficiency). 
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Figure  3.  Simulation  MFM  Link 


Table  2.  SUMMARY  OF  CODES  IN  SIMULATION 


Code  (n,k) 


Block  Length  (n) 


Information  Symbols 
(k) 


Parity  Symbols  (n  -  k) 


cU 


(15,12) 


15 


»/(Bits;  Hz)  3.0 


*  t  Errn  ,  ,•  SEC/DED 

Correction/ Detection  ' 


*  can  correct  some  triple  burst  errors 


(15,10) 


15 


(31.2S) 


31 


B.  CHANNEL  CODING  AND  DECODING 

The  16-QAM  constellation  consists  of  two  magnitude  rings  and  eight  phases. 
The  magnitude  rings  are  of  amplitudes  1  A'  and  2 A',  with  the  value  of  X  variable. 
Data  from  simulation  runs  was  collected  with  X  =  100.  The  phases  employed 
were  offset  ±22.5*  from  the  axes.  Figure  4  is  a  diagram  of  the  16-QAM 


20 


constellation  and  Appendix  A  demonstrates  the  one-to-one  relation  between 
Galois  field  elements  and  QAM  constellation  elements  for  both  16-QAM  and 
32-QAM.  Decoding  for  this  constellation  is  accomplished  for  each  sample  in  the 

information  band  of  the  FFT  as  follows: 

1.  Determine  the  quadrant  of  the  received  complex  signal  by  the  sign  of  its  real 
and  imaginary  parts. 

2.  Determine  which  of  the  two  phases  in  the  quadrant  is  closest  to  the  signal. 
This  is  done  by  comparing  the  absolute  value  of  the  real  pau  to  that  of  the 
imaginary  part.  If  that  of  the  real  part  is  greater,  choose  the  phase  closest  to 
the  real  axis.  If  not,  choose  the  phase  closest  10  the  imaginary  axis. 

3.  Finally,  estimate  the  correct  magnitude:  compare  the  received  magnitude 
with  1.5A'.  If  the  received  magnitude  is  greater  than  1 .5 A’,  assume  that  the 
received  signal  is  actualiy  the  constellation  element  on  the  outer  ring.  If  not, 
assume  it  is  the  dement  on  the  inner  ring. 

Two  32-QAM  constellations  were  used  in  the  simulation.  The  first  is  Gray 
coded  and  consists  of  four  magnitude  rings  and  eight  phases.  The  magnitude 
rings  are  of  amplitudes  LA,  2X.3X  and  4A\  once  again  with  X  variable.  Data 
for  32-QAM  was  collected  with  X  =  20.  The  eight  phases  are  the  same  as  for 
16-QAM,  ±  22.5’  offset  from  the  axes.  Figure  5  is  a  diagram  of  this  Gray-coded 
constellation.  The  second  32-QAM  constellation  is  identical  to  the  first  except 
that  the  rings  of  magnitude  2A'  and  4A'  are  shifted  clockwise  22.5”  to  increase  the 
minimum  distance  between  constellation  elements.  Figure  6  is  a  diagram  of  the 
modified  32-QAM  constellation.  Channel  decoding  for  both  32-QAM 
constellations  uses  a  different  algorithm  than  16-QAM,  one  which  finds  that  ele¬ 
ment  of  the  constellation  that  is  closest  in  distance  to  the  received  signal  (a  form 
of  maximum-likelihood  decoding).  It  is  described  as  follows,  with  rx  being  a  re¬ 
ceived  complex  sample  in  the  information  band  of  the  FFT  and  q, ,  i  =  0,1,. ..,31 
being  the  elements  of  the  32-QAM  constellation: 
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figure  4.  16-QANI  Constellation 


1.  Let  ymm  =  J(Rc{rx)  -  Re(q0))2  +  (Imfo)  -  Im(<70))2  ,  the  length  of  a  line  seg¬ 
ment  drawn  between  rx  and  q0  on  the  real-imaginary  axis. 

2.  Compare  ymn  to  y,  =  J{R e(rx)  -  Re{q,))2  +  (Im(rx)  -  \m{q,))2  for  each  suc¬ 
cessive  element  of  the  constellation.  If  y,  is  less  than  ymin,  replace  the  value 
of  )'ran  with  the  value  of  yn  store  the  identity  of  the  constellation  element 
which  gave  the  smaller  value  (denoted  by  the  subscript  i ),  and  try  the  next 
element. 

3.  After  comparison  to  all  elements  of  the  constellation,  i  identifies  the  closest 
element.  Assume  that  the  element  i  is  the  actual  symbol. 


The  QPSK  constellation  consists  of  four  phases:  one  each  at  ±  45’  and 
±  135'.  Data  was  collected  with  these  points  having  magnitude  ten.  Channel  de¬ 
coding  consists  of  determining  the  quadrant  of  the  received  samples  based  on  the 
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Imag. 


Figure  5.  Gray  Coded  32-QAM  Constellation 

signs  of  their  real  and  imaginary  parts.  The  received  phase  is  assumed  to  be  that 
which  lies  within  the  received  quadrant. 

C.  SUPPORTING  FILES  AND  PROGRAMS 

Prior  to  writing  an  MFM  link  program,  two  data  files  were  created  ("GF16 
DATA"  and  "GF32  DATA")  which  contain  the  binary  representations  of  the  el¬ 
ements  of  GF(2“)  and  GF(25),  respectively.  These  files  form  the  basis  of  the 
Reed-Solomon  computations  in  the  link.  In  addition,  three  data  files  ("QAM  16", 
"QAM32"  and  "QAM32A")  were  created  which  outlined  the  QAM  constellations 
used,  containing  magnitude  and  phase  information  on  each  element. 
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Imag. 


Figure  6.  Modified  32-QAM  Constellation 

Several  other  programs  were  written  to  support  the  operation  of  MFMLINK. 
The  program  "GALOIS"  reads  the  w-bit  representations  of  a  Galois  field  and 
outputs  the  GF(2m)  addition  table.  GALOIS  adds  each  element  bit-by-bit 
(modulo-2)  and  then  compares  the  result  against  each  element  of  the  input  data 
file  to  determine  which  element  of  GF(2m)  it  corresponds  to.  This  addition  table 
is  output  to  the  data  file  "GFA16"  or  "GFA32"  as  appropriate.  The  program 
"BITDIFF"  reads  in  the  w-bit  representations  of  GF(2m)  and  compares  each  ele¬ 
ment  bit-by-bit  to  determine  the  number  of  places  in  which  each  differs  for  use 
in  calculating  the  number  of  bit  errors.  The  output  of  BITDIFF  is  a  GF(2m) 
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"difference"  table,  either  "GFD16"  or  "GFD32".  The  16-elemcnt  difference  table 
for  the  QPSK  symbols  was  entered  separately  in  the  file  "QPSKD1".  "TRACE" 
reads  in  the  Galois  field  addition  tables  created  by  "GALOIS"  and  computes  the 
trace  of  each  element  (sec  eqn.  20).  Finally,  "QAMCOM"  reads  in  the  magnitude 
and  phase  information  of  each  QAM  constellation  and  calculates  the  corre¬ 
sponding  real  and  imaginary  components.  This  output  is  placed  in  data  file 
"QAMRI16",  "QAMRI32  or  "QAMRI32A",  for  use  in  channel  encoding  and 
decoding.  The  real  and  imaginary  components  for  the  QPSK  constellation  were 
entered  separately  in  the  data  file  "QPSKRI". 

The  outputs  of  MFMLINK  are  written  into  the  data  files  "RSRSLTS"  and 
another  file  with  a  variable  name,  depending  on  the  coding  scheme  used. 
RSRSLTS  contains  information  on  code  chosen,  total  symbol  and  bit  errors, 
number  of  bauds  run,  and  (if  desired)  the  actual  random  data  symbols  input, 
transmitted  code  words,  received  signals,  and  corrected  code  words.  The  file  with 
the  variable  name  contains  only  the  SNRs  used  in  the  run  and  resulting  bit  error 
rates  for  use  in  plotting,  and  is  named  according  to  the  code  and  baud  type  of  the 
run  (e.g.,  "C0BT1"  corresponds  to  code  "0"  (uncoded  16-QAM)  and  baud  type 
"1"). 

D.  MFMLINK 

MFMLINK  is  subroutine  based,  for  e„se  in  altering  parameters  or  algo¬ 
rithms.  One  of  the  difficulties  in  adapting  a  finite  Galois  field  into  a  software 
simulation  is  the  use  of  non-binary  field  arithmetic,  the  addition  and  multipli¬ 
cation  of  the  powers  of  primitive  element  a.  In  the  program  (and  supporting  files), 
this  was  handled  by  identifying  field  elements  in  terms  of  its  power  of  a  plus  one, 
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except  for  the  zero  element  which  is  always  "0"  (i.e.,  0  -+  0,  1  =  a0  -*  1,  a1  -»  2, 
etc.).  With  this  method  of  element  identification  for  GF(2m),  field 
multiplication  division  of  elements  other  than  zero  is  carried  out  by 
addition,  subtraction  modulo  n  =  2m  -  1. 

A  program  run  starts  with  the  user  choosing  the  type  of  coding  scheme  and 
baud,  and,  if  error  control  is  utilized,  whether  code  words  with  more  errors  de¬ 
tected  than  the  code  can  correct  should  be  erased.  SNRs  and  the  number  of 
bauds  for  each  SNR  are  fixed  in  the  program  and  must  be  edited  in  order  to 
change.  Once  the  user  inputs  these  parameters,  subroutine  "BAUDS"  reads  the 
FFT  size  and  its  logarithm  base  two,  ku  k2,  and  the  number  of  code  words  per 
baud  (when  error  control  is  used).  Subroutine  "DATA"  then  reads  in  the  Galois 
addition  and  difference  tables,  element  traces,  and  QAM  constellation  charac¬ 
teristics.  Once  this  background  data  is  read,  the  run  commences. 

If  an  uncoded  case  is  chosen,  subroutine  "NOCODE"  provides  random 
Galois  field  elements  as  data  and  the  corresponding  QAM  or  QPSK  constellation 
elements,  and  loads  the  inverse-FFT  array  for  transmission.  For  coded  cases,  the 
main  program  determines  the  random  Galois  field  data  for  a  code  word,  which 
is  then  encoded  in  Reed-Solomon  systematic  form  (with  subroutine  "ENCODE"). 
This  process  is  repeated  until  the  baud  is  filled  with  Reed-Solomon  code  words, 
with  each  tone  representing  one  element  of  the  Galois  field  in  use.  The  main 
program  then  loads  the  corresponding  inverse-FFT  array  for  transmission. 

Subroutine  "BITREV"  places  the  array  in  bit-reversed  order  for  use  in  the 
subroutine  "1NVFFT"  which  takes  tne  inverse-FFT  of  the  transmission  array 
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using  the  decimation-in-time  algorithm.  Both  BITREV  and  1NVFFT  were  copied 
from  the  digital  signal  processing  files  of  Professor  Robert  D.  Strum  at  NPS. 

Once  a  baud  has  been  transmitted,  noise  is  added  to  each  sample  of  the  time 
domain  signal.  First,  signal  power  is  calculated  from  the  sum  of  the  squares  of  the 
time  domain  signal  samples  (a  real-valued  sequence).  Then  noise  power  (or  var¬ 
iance,  since  the  noise  is  fixed  to  have  zero  mean)  is  calculated  from 

A-,-1 

,  I  Tx(if 

2  ___  1  i — o 

C7W_  ('Yr —  1 )  1q[SSR-i,)I\0.Q  1  ’ 

where  cl  is  the  noise  power,  Ns  is  the  FFT  size,  Tx  is  the  inverse- F FT  array 
(transmitted  signal),  and  ??  is  given  by 


(30) 


The  quantity  >/  corrects  the  value  of  SNR  to  be  the  ratio  of  signal  power  to  noise 
power  in  the  frequency  band  from  /c,  to  /c2,  rather  than  having  the  signal  power 
concentrated  in  this  frequency  interval  and  noise  power  spread  uniformly 
throughout  the  spectrum.  Once  the  noise  variance  is  calculated,  subroutine 
"GAUSS"  calculates  random  samples  of  AWGN  with  zero  mean  and  variance 
cl  to  add  to  each  sample  of  the  time  domain  signal,  forming  the  received  signal 
samples. 

The  subroutines  BITREV  and  "FFT"  (FFT  is  also  copied  from  Professor 
Strum's  files)  function  as  the  MFM  receiver.  Each  received  complex  sample  in  the 


27 


information  band  from  /c,  to  k2  is  then  channel  decoded  into  an  element  of  the 
appropriate  Galois  field  by  the  subroutines  "WHAT16"  (for  16-QAM), 
"WHAT32"  (for  32-QAM)  or  "WHATQP"  (for  QPSK).  For  uncoded  cases,  these 
decoded  symbols  are  compared  with  the  input  symbols  to  determine  symbol  and 
bit  errors.  For  runs  with  error  control,  after  channel  decoding  each  estimated 
code  word  is  then  decoded  using  the  appropriate  algorithm  of  Chapter  111.  First, 
each  estimated  codeword's  syndrome  is  calculated  (with  subroutine  "SYNDRO"). 
If  any  errors  are  detected,  either  subroutine  "DECONE"  (for  SEC  codes)  or 
"DECTWO"  (for  DEC  codes)  are  called  to  correct  errors,  if  possible.  If  code  word 
erasure  was  called  for  by  the  user,  these  decoding  subroutines  call  "RUBOUT" 
to  erase  code  words  that  cannot  be  corrected.  After  Reed-Solomon  decoding,  the 
input  code  words  are  compared  to  the  corrected  output  to  determine  the  number 
of  symbol  and  bit  errors.  Finally,  subroutine  "OUTPUT"  sends  pertinent  output 
information  to  the  data  file  RSRSLTS. 

E.  RESULTS 

The  primary  results  of  the  simulation  are  displayed  in  the  plots  of  bit  error 
rate  versus  SNR  and  bit  error  rate  versus  EB/N0  .  At  lower  values  of  SNR  (or 
Eb/iX0),  accurate  results  were  obtained  by  performing  runs  with  10,000-20,000 
bits.  At  higher  values  of  SNR,  as  the  number  of  errors  diminished,  over  100,000 
bits  were  required.  The  results  of  the  simulation  runs  are  plotted  in  the  eight 
graphs  included  in  the  following  pages.  The  graphs  are  in  pairs  and  reflect  the 
bit  error  rates  as  a  function  of  either  SNR  or  EB/N0  .  The  plots  with  Eb/iX0  are 
derived  from  those  with  SNR  based  on  the  following  relationship: 
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where  SNRus  is  the  wide  band  SNR  (that  used  as  an  input  to  the  program)  and 
rn  is  the  number  of  bits  per  tone  for  a  particular  modulation  scheme.  The  calcu¬ 
lation  of  EbI A’o  includes  the  energy  in  all  bits  carried  on  a  tone,  both  information 
and  parity,  since  the  bit  error  calculations  include  errors  in  both  data  and  parity 
symbols.  Figure  7  is  a  graph  of  the  uncoded  modulation  schemes'  performances 
on  the  link  in  terms  of  Bit  Error  Rate  vs.  SNR  and  also  shows  the  same  per¬ 
formance  in  terms  of  EBI A"0  .  Figure  S  contains  similar  plots  showing  the  per¬ 
formance  of  the  two  Recd-Solomon  codes  used  with  16-QAM  against  that  of 
uncodcd  16-QAM  and  uncoded  QPSK.  Figures  9  and  10  show  the  performance 
of  the  two  Rced-Solomon  codes  used  with  both  32-QAM  constellations  against 
their  respective  uncoded  cases.  A  summary  of  the  BER  performances  of  the  error 
control  codes  over  the  corresponding  uncoded  cases  is  contained  in  Table  3  fol¬ 


lowing  the  performance  graphs. 


Figure  7.  Uiicoded  Bit  Error  Rates:  Top:  HER  vs.  SNR.  Bottom:  BCR  vs.  ZT,//V0 
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Figure  9.  Gray  Coded  32-QAM  Bit  Error  Rates:  Top:  BER  vs.  SNR.  Bottom: 
BER  vs.  EbI iV0 . 


32 


Figure  10.  Modified  32-QAM  Bit  Error  Rates:  Top:  BLR  vs.  SNR.  Bottom:  BF.R 
vs.  r„i A”0 . 
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Table  3.  ERROR  CONTROL  PERFORMANCE  SUMMARY 


16-QAM 

Code 

(15.12) 

(15.10) 

BER 

IO'3 

10-4 

BB 

io-3 

io-4 

io-5 

Coding 

Gain 

(approx.) 

2.5 

2.5 

2.5 

4.0 

4.7 

5.0 

Gray  Coded  32-QAM 

Code 

(31.28) 

(31.26) 

BER 

io-3 

io-4 

io-5 

io-3 

10-4 

Coding 

Gain 

(approx.) 

1.4 

1.5 

1.5 

2.4 

2.6 

2.6 

Modified  32-QAM 

Code 

(31.28) 

(31.26) 

BER 

10  3 

io-4 

io-5 

io-3 

io-4 

10-5 

Coding 

Gain 

(approx.) 

2.1 

2.1 

2.5 

2.9 

3.5 

3.8 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  simulation  MFMLINK  is  the  First  implementation  of  error  control  cod¬ 
ing  in  multi-frequency  modulation  at  NPS.  As  such,  the  codes  utilized  may  not 
be  the  optimal  codes  for  the  given  link  parameters,  but  do  provide  good  insight 
into  the  benefits  of  error  control  coding  for  M  KM  - 

First,  the  enhanced  error  detection  capabilities  of  the  modified  Reed-Solomon 
codes  in  the  simulation  support  the  introduction  of  some  type  of  Automatic  Re¬ 
peat  Request  (ARQ)  strategy  into  the  MFM  environment.  More  specifically,  a 
type  I  hybrid  ARQ  scheme  is  called  for  as  described  in  Reference  2.  Using  this 
scheme  with  MFM,  any  correctable  errors  which  are  detected  are  corrected  and 
if  uncorrectable  errors  are  detected,  the  receiver  requests  retransmission  of  par¬ 
ticular  erroneous  codewords  or  bauds.  The  proper  combination  of  ARQ  and  FEC 
can  provide  a  higher  information  throughput  than  one  with  ARQ  alone,  and 
provide  higher  system  reliability  than  one  with  FEC  alone  [Ref.  2:  p.478]. 

Second,  although  the  decoding  time  for  any  given  codeword  or  baud  was  not 
measured,  the  use  of  the  decoding  algorithms  presented  for  the  modified  Reed- 
Solomon  codes  (which  are  solely  based  on  the  calculated  syndromes)  are  simpler 
to  implement  in  software  than  the  iterative  algorithms  for  conventional  codes  and 
felt  to  be  less  time  consuming  to  run.  The  loss  in  information  rate  due  to  surren¬ 
dering  an  information  symbol  to  parity  is  well  worth  the  time  saved  in  the  de¬ 
coding  process. 
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Finally,  based  on  the  graphs  presented  at  the  end  of  Chapter  4,  the  coding 
gains  obtained  at  bit  error  rates  of  10~3  and  10-4  for  each  code  has  proven  the 
value  of  the  codes  used.  Large  gains  were  evident  using  error  control  with 
16-QAM  modulation  at  a  moderate  cost  in  code  rate.  Smaller  gains  were  evident 
using  error  control  with  32-QAM  modulation,  particularly  with  the  Gray-coded 
constellation,  but  at  only  a  slight  reduction  in  code  rate.  In  any  case,  the  loss  in 
code  rate  due  to  the  addition  of  parity  symbols  for  error  control  is  worth  the 
coding  gains  achieved. 

Since  this  simulation  was  the  first  venture  into  MFM  error  control  imple¬ 
mentation  at  NPS,  several  alterations  to  the  simulation  may  provide  improved 
error  control  performance.  One  might  be  to  change  the  methods  of  channel  de¬ 
coding.  particularly  in  the  case  of  32-QAM.  The  approach  used  with  16-QAM 
may  have  improved  performance  over  the  maximum-likelihood  approach. 
Constellation  geometries  with  elements  farther  apart  may  improve  performance, 
as  seen  in  the  results  obtained  with  the  32-QAM  constellations.  Another  alter¬ 
ation  might  be  to  vary  the  modulation  technique  used  with  the  error  control  codes 
in  the  simulation.  For  example,  using  QPSK  and  the  length  15  Reed-Solomon 
codes  with  each  tone  carrying  two  bits,  two  consecutive  tones  would  represent  one 
element  of  GF(24)  and  30  tones  would  represent  a  complete  codeword.  Interleav¬ 
ing  the  codeword  symbols  between  bauds  is  another  possible  alteration,  e.g., 
having  a  16  tone  information  band  with  each  tone  representing  one  symbol  from 
16  separate  codewords  would  correspond  to  the  transmission  of  16  codewords  in 
a  15  baud  packet.  Other  methods  of  error  control  are  also  of  interest,  particularly 
the  use  of  convolutional/trellis  codes  and  concatenated  codes. 
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APPENDIX  A.  BINARY  REPRESENTATIONS  OF  GALOIS 


TabJe  4.  ELEMENTS  OF  GF(  2* ) 


Primitive  Polynomial: 

P{x)  =  1  +  x  +  x4 

Element 

Binary 

Representation 

0 

0000 

1000 


0100 


0010 


0001 


1100 


01 10 


0011 


!  101 


1010 


0101 


1110 


0111 


mi 


1011 


1001 


Table  5.  ELEMENTS  OF  GF(  25 ) 


Primitive  Polynomial: 

P(x)  =  \  +  X2  +  X5 

Element 

Binary 

Representation 

o 

00000 
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APPENDIX  B.  SUPPORTING  PROGRAMS 


A.  GALOIS 


c 

c 

c 

c 

c 

c 


Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  ‘V  Vr  -V  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  *  “V  Vr  Vr  *  *  Vr  Vr  -jV  Vr  Vr  *  Vr  Vr  Vr  Vr  Vr  Vr  Vr  -sV  Vr  Vr  Vr 


THIS  PROGRAM  WILL  READ  THE  M-BIT  REPRESENTATIONS  OF  GF(2**M)  FROM  * 
A  DATA  FILE,  THEN  CALCULATE  THE  ADDITION  TABLE  ( (N+1)*(N+1) )  AND  * 
WRITE  THE  RESULTS  INTO  A  DATA  FILE  GFA(N+1).  * 

* 


*  Vc  Vr  Vr  Vr  'V  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  ■ 


■  Vc  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  *  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  *  Vr  Vr  *  *  Vr  Vr  Vr 


INTEGER  P ( 0 :  15),G(0:  15,0: 7),A(0: 15,0: 15,0: 7),GFA(0: 15,0: 15) 
INTEGER  S  , N , M 

CALL  EXCMS( ' FILEDEF  8  DISK  GF16  DATA  (PERM') 

OPEN( 7, FILES' GFA16' ,ACCESS=' SEQUENTIAL' , FORM= ' FORMATTED ' ,STATUS= 
+  'UNKNOWN') 

OPEN( 8 ,FILE=' GF16 ' ,ACCESS=' SEQUENTIAL' ,FORM=' FORMATTED' , STATUS 
+  =’ UNKNOWN') 

N=15 
M=4 

C  READ  IN  THE  M-BIT  REPRESENTATIONS  FROM  DATA  FILE 


DO  10  I  =  0 ,N 

READ( 8,99)  P( I) , (G( I , J) , J=0 ,M-1) 
WRITE(5 ,99)  P( I) , (G( I , J) ,J=0,M-1) 
10  CONTINUE 

99  FORMATf IX, 12,412) 

C  ADD  THE  M-BIT  ELEMENTS.  .  . 

DO  20  1=0, N 

DO  30  J=0,N 

DO  40  K=0 ,M- 1 

S=G( I ,K)+G( J,K) 

A( I , J ,K)=MOD( S , 2) 

40  CONTINUE 

30  CONTINUE 

20  CONTINUE 


C  NOW  FIGURE  OUT  WHICH  POWER  OF  ALPHA  THESE  M  BITS  REPRESENT. .  . 

DO  50  J=0,N 

DO  60  J=0,N 

DO  70  K=0 ,N 

DO  80  L=0 ,M-1 

IF  (  A(  I ,  J ,  L) .  NE.  (  G(  K ,  L)  )  )  THEN 
GO  TO  70 
END  IF 


oooooooonnoooooo 


80  CONTINUE 

GFA( I , J)=P(K) 

GO  TO  60 

70  CONTINUE 

60  CONTINUE 

50  CONTINUE 

C  NOW  WRITE  THESE  POWERS  OF  ALPHA  INTO  THE  ADDITION  TABLE.  .  . 


DO  90  1=0, N 

WRITE( 7 ,9001  ( GFA( I , J) , J=0 , 15) 

THE  FOLLOWING  LINES  ARE  FOR  LARGER  FIELD  CALCULATIONS.  .  . 
WRITE( 7 , 900)  (GFA(I,J),J=16,31) 

WRITE (7 ,900)  (GFA(I,J),J=32,47) 

WRITE( 7 , 900)  (GFA( I , J) , J=48 ,63) 

WRITEf  7 ,900)  (GFA(I.J) ,J=64,79) 

WRITE( 7 , 900)  ( GFA( I , J) , J=80 , 95 ) 

WRITE( 7,900)  (GFA(I , J) , J=96, 111) 

WRITE( 7,900)  (GFA(I,J),J=112,127) 

WRITE ( 7,900)  (GFA(I.J), J=128,143) 

WRITEC 7,900)  (GFA(I.J) ,J=144,159) 

WRITE( 7 ,900)  (GFA( I , J)  ,J=160,175) 

WRITE( 7, 900)  (GFA(I, J) ,J=176,191) 

WRITEC 7, 900)  (GFACI, J) ,J=192,207) 

WRITEC 7, 900)  (GFA(I, J),J=208, 223) 

WRITEC 7,900)  (GFA(I, J),J=224, 239) 

WRITEC  7,900)  ( GFA( I , J ) , J=240 , 255 ) 

90  CONTINUE 

900  FORMATC 16( IX, 12) ) 

ENDFILEC  7) 

CLOSEC  7) 

END 
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B. 

c* 

c 

c 

c 

c* 


BITD1FF 


-,V  V?  V?  «'r  “V  Vr  iV  vV  ~V  Vc  Vr  Vc  Vr  *V  •&  Vr  “V  Vr  V.*  V c  ☆  *  “V  Vr  sV  Vr  “V  %■  *' 


•  • 


•** 


THIS  PROGRAM  READS  IN  THE  BIT  REPRESENTATION  OF  A  GALOIS  FIELD  AND  * 
CALCULATES  THE  #  OF  BITS  IN  WHICH  EACH  ELEMENT  DIFFERS.  THEN  IT  * 
WRITES  THE  RESULTS  IN  A  TABLE  FOR  USE  IN  FIGURING  BIT  ERROR  RATE  * 


INTEGER  ALPHA ( 0:  31),BIT(0:  31 , 0:  5)  ,GFD( 0:  31,0: 31) ,N,ROWS,DIFF 

CALL  EXCMS( 'FILEDEF  9  DISK  GF32  DATA  (PERM') 

OPEN  ( 7 ,FILE=' GFD32 ' , FORM= ' FORMATTED ' ,ACCESS=' SEQUENTIAL' , 

+  STATU S= ' UNKNOWN ' ) 


N=31 

P0WER=5 

R0WS=(N+1)/16 

DO  10  1=0, N 

READ(9,*)ALPHA(I) ,(BIT(I,J) ,J=0, POWER-1) 
10  CONTINUE 

DO  20  1=0, N 
DO  30  J=0,N 
DIFF=0 

DO  40  K=0, POWER- 1 

IF  ( BIT( I ,K) . NE. BIT( J ,K) )  THEN 
DIFF=DIFF+1 
END  IF 

40  CONTINUE 

GFDf I , J)=DIFF 
30  CONTINUE 

20  CONTINUE 

DO  50  1=0, N 

WRITE ( 7 , 1000) (GFD( I , J) , J=0 , 15) 

WRITE( 7 , 1000)(GFD( I , J) ,J=16,31) 

50  CONTINUE 

1000  FORMATC 16(12, IX)) 

END 
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on  o  o  o  o  o 


C.  TRACE 


»*«•  <.'•  »*«  »*«  Vr  Vr  Vc  Vr  V*  V%  #V  *Wr  «Wr  Vr  *  V  V'V*  *V  *  V  #V  Vw « V%  Vr  V«  V«  V*  #W*  Vc 

THIS  IS  THE  TRACE  FINDER  PROGRAM  * 


'f***«fif********Vf*# 


\ 

INTEGER  A,B,C,D,GFA(0: 31,0: 31), N 

CALL  EXCMS( ' FILEDEF  9  DISK  GFA32  FILE  (PERM') 

OPEN( 6 , FILE= ' TR32 ' , FORM® ' FORMATTED ' ,ACCESS=' SEQUENTIAL’ , 
+  STATUS® ' UNKNOWN ' ) 


READ  IN  GF(2**M)  ADDITION  TABLE... 
N=31 

DO  10  1=0,31 

READ( 9 ,*)  (GFA( I , J) , J=0 , 15 ) 
READ( 9 ,*)  (GFA( I , J) , J=16 , 31) 
10  CONTINUE 

C . 

VRITEC6.1001) 

1001  FORM AT ( ' 0 ' , IX, ' 0 ’ ) 

DO  55  1=2,31 

A=M0D( 2*( I - 1) , N ) + 1 
A=GFA( I ,A) 

B=M0D(4*( 1-1) ,N)+1 
B=GFA( A,B) 

C=M0D(  8*(  I  - 1)  ,N)+1 
C=GFA( B , C) 

D=M0D( 16*(I-1) , N ) + 1 
WRITE( 6 , 1002)GFA(C,D) 

55  CONTINUE 

1002  FORMATC IX, 12) 

END 
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D.  QAMCOM 


QAM  COMPUTER-  THIS  PROGRAM  READS  IN  MAGNITUDE  AND  PHASE  INFO 
ON  A  QAM  CONSTELLATION  AND  CONVERTS  IT  TO  A  REAL  AND 
IMAGINARY  PART  WHICH  IS  STORED  IN  A  DATA  FILE  "QAMRl" 
WHICH  IS  READ  BY  THE  MFM  LINK  PROGRAM. 


INTEGER  ALPHA( 0: 31),MAG(0: 31) ,N 

REAL  RE( 0:  31),IM(0: 31) ,PI ,PHASE(0:  31) 

CALL  EXCMS( 1 FILEDEF  9  DISK  QAM32A  DATA  (PERM') 

OPEN( 6 , FI LE= ' QAMRl 32A ' , FORM= ’ FORMATTED ' , ACCESS= ' SEQUENTIAL ' 
+  STATUS*' UNKNOWN’) 


N=31 

PI =4.  0*ATAN( 1.0) 


DO  1  1=0, N 

READ(9 ,*)ALPHA( I) ,MAG(  I) ,PHASE(  I) 
PHASEC I )=PHASE( I )*PI/ 180.  0 
RE ( I )  =MAG(  I )  •'•'COS ( PHASE ( I ) ) 

IM(  I )=MAG( I )*SIN( PHASE( I ) ) 

WRITE ( 6 ,*)RE( I ) , IM( I ) ,MAG( I )/l. 0 
1  CONTINUE 

ENDFILEC6) 

CLOSE( 6) 

END 


43 


ooooooooooonoooooooooooonooooooo 


APPENDIX  C.  SIMULATION  PROGRAM:  MFMLINK 


V?  Vr  Vc  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vc  'V  Vr  Vc  ‘V  Vc  Vr  ‘.V  Vr  Vr  Vr  Vr  Vc  Vc  Vr  Vc  Vc  Vr  Vr  Vc  Vc  Vr  Vc  Vc  ‘V  Vr  Vr  Vr  Vc  ‘V  Vr  Vr  Vr  Vr  ‘V  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vc  Vc  Vr  Vc  Vr  Vr  Vr  Vr  Vr  Vr  "V  Vr  Vr  Vc  ■jV  Vr  -.V 
*/fVf*Vf*VrVT*Vr*Vc****VriV*****Vf*Vf*^VVtVrVcV?^VVtVrVfVoV^VVrVrVc**Vr^V^V*Vc**>V*VfVf“V***-}V**Vr,3V**^V’jV->V’;V’jV 
Vc  VcVr 

*  MFM  COMMUNICATION  LINK 


-THIS  SIMULATES  AN  MFM  COMMUNICATION  LINK  WHICH  IS 
CORRUPTED  BY  AWGN.  IN  DEALING  WITH  GALOIS  FIELDS 
IT  WAS  NECESSARY  TO  TREAT  THE  ELEMENTS  IN  TERMS 
OF  THEIR  POWER  OF  ALPHA  PLUS  ONE,  EXCEPT  FOR  THE 
ELEMENT  "0"  WHICH  IS  ALWAYS  WRITTEN  AS  "0". 

Vr  Vr  Vc  Vc  Vc  Vr  V.-  Vc  Vc  Vc  Vc  Vr  'V  *V  Vr  Vr  Vr  Vr  Vr  Vr  Vc  *  V  Vc  Vc  Vr  Vr  Vc  Vr  Vr  -.V  Vr  Vr  Vc  Vr  Vr  V?  Vc  Vr  Vc  *  ‘V  Vr  ‘V  Vc  Vr  Vr  Vr  Vr  Vr  Vr  Vc  Vr  Vr  Vc  Vr  Vr  Vr  Vr  Vr  Vc  Vr  -V  Vr  Vr  ■> 

-THE  IMPORTANT  ARRAYS /VARIABLES  ARE  AS  FOLLOWS: 

ARRAYS:  D . RANDOM  DATA  VECTOR  (CODED  CASES) 

DIN . RANDOM  DATA  VECTOR  (UNCODED  CASES) 

CW . REED -SOLOMON  CODEWORD  FOR  D 

TX . ANALOG  TRANSMISSION  (OUTPUT  OF  IN'V-FFT) 

AWGN . NOISE  VECTOR 

RX . ANALOG  RECEPTION  (TX  +  AWGN,  OUT  OF  FFT) 

R . RECEIVED  VECTOR 

DOUT . RECEIVED  VECTOR  (UNCODED  CASES) 

E . ESTIMATED  ERROR  IN  R 

V . ESTIMATED  CODEWORD(R  +  E) 

DHAT . ESTIMATED  DATA  FROM  R 

GFA( I, J). ...  GALOIS  FIELD  SUM  OF  SYMBOLS 

GFD( I, J). ...  NUMBER  OF  BITS  THAT  I  AND  J  DIFFER  IN 

S . RECEIVED  VECTOR  SYNDROME 


VrVrVrVrVrVrVrVrVrVcVcVrVrVrVr  VrVrVrVrVr*Vr*VrVrVrVr?VVrVoVVrVrVrVrVrVrVrVrVr*VrVrVrVrVrVrVrVr**VrVrVr^VVrVr^WrVrVr**VrVrVrVrVrVrVr 


INTEGER  B , C , K , BYTERR , ERRORS , FIXED, BDS, RUB , NLESSK , MAX 

INTEGER  N,M, ITER, ERASE, DATERR,CCWW( 17,0:  31) ,DIN(0:  255) ,DOUT(0:  255) 

INTEGER  S( -2:  2),GFA(0: 31,0: 31),R(0: 31),V(0: 31),E(0: 31) 

INTEGER  TRACE ( 0: 31) ,T4(0: 31) ,D(0: 31) ,CW(0: 31) ,DD( 17,0:  31) 

INTEGER  GAMMA 1 , GAMMA 2 , GAMMA 3 , ZEROS ,TR I ERR, DUOERR ,UNOERR , OOPS 
INTEGER  DHAT( 0:  31) ,BDTYPE ,CODE ,K1 ,K2 , PTS ,CPB , PTR, BAUDNO 
INTEGER  GFD32( 0: 31,0: 31) ,DBERR,CWBERR,GFD16( 0: 15,0: 15) 

INTEGER  QPDIFF(0: 3,0: 3) 

REAL  Dl( 0:  31) ,D0(0: 255) ,QAMRE( 0: 31) ,QAMIM(0: 31) , QAMMAG( 0: 31) 

REAL  MEAN , NVAR , AWGN ( 0 : 4095 ) ,TEMP1( 0:  4095 ) ,TEMP2( 0:  4095) 

COMPLEX  QAMCOM( 0: 31) ,TX(0: 4095) ,RX(0: 4095) ,T(0: 4095) 

REAL  SUM , SUMSQ , SNR , NF IX , QPRE ( 0 : 5) ,QPIM(0: 3) 

CALL  EXCMS( ’FILEDEF  9  DISK  GFA16  DATA  (PERM') 

CALL  EXCMS( 'FILEDEF  7  DISK  TR16  DATA  (PERM)') 

CALL  EXCMS( 'FILEDEF  8  DISK  QAMRI16  TILE  (PERM') 

CALL  EXCMS( 'FILEDEF  1  DISK  GFD16  FILE  (PERM') 

CALL  EXCMS( 'FILEDEF  11  DISK  GFA32  FILE  (PERM’) 
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CALL  EXCMS( ' FILEDEF  12  DISK  TR32  FILE  (PERM)') 
CALL  EXCMS( 'FILEDEF  13  DISK  Q AMR 1 32 A  FILE  (PERM') 
CALL  EXCMS( 'FILEDEF  31  DISK  GFD32  FILE  (PERM') 
CALL  EXCMS( 'FILEDEF  14  DISK  QPSKRI  FILE  (PERM') 
CALL  EXCMS( 'FILEDEF  34  DISK  QPSKDI  FILE  (PERM') 


OPEN’(  6  ,FILE='  RSRSLTS  '  , FORM* ' FORMATTED ’  .ACCESS*' SEQUENTIAL'  , 

+  STATUS*' UNKNOWN’ ) 

OPEN( 33 ,FILE='  C2C1NEW' , FORM* ' FORMATTED ' , ACCESS*' SEQUENTIAL' , 
+  STATUS* ' UNKNOWN ' ) 


C 

c 


SELECT  WHICH  CODE  DESIRED... 
WRITE (5.*)'  ' 


413 


WRITE( 5 ,*) ' 1. 
WRITE( 5 ,*) : 
WRITE( 5 ,*) ' 
WRITE( 5 ,*) ' 

WR I TE ( 5 , * ) ' 
WRITE (5,*) ’ 
WRITE  (5,  •••-)' 
WRITEC  5,  ••'■•)’ 
READ(5  ,';r)CODE 


SELECT  CODE  TYPE.  .  . 

0  =  UNCODED  (16 -QAM)' 

1  =  (15,10)  REED-SOLOMON  CODE 

2  =  (31,26)  REED  SOLOMON  CODE 

3  =  (15,12)  REED  SOLOMON  CODE 

4  =  (31,28)  REED  SOLOMON  CODE 

5  =  UNCODED  (32 -QAM)' 

6  =  UNCODED  QPSK' 


( 16-QAM) ' 
( 3 2 -QAM) ' 
( 16-QAM) ' 
(32 -QAM)’ 


IF  ( ( CODE.  LT.  0) .  OR.  ( CODE.  GT. 6) )  THEN 

WR I TE ( 5 , * ) ' BAD  CHOICE  OF  CODE --TRY  AGAIN?' 

GO  TO  413 
END  IF 

IF  (CODE. EQ. 6)  THEN 
WRITE(  5  ,•'•')  ’  ’ 

21  WRITE(  5  ,■'•■) '  2.  SELECT  BAUD  TYPE  (1-5)’ 

READ( 5 ,*)3DTYPE 

IF  ( ( BDTYPE.  GT.  5 ) .  OR.  ( BDTYPE.  LT.  1 ) )  THEN 
GO  TO  21 
END  IF 
END  IF 

IF  ((CODE.  EQ.  1).  OR.  (CODE.  EQ.  0).OR.  (CODE.  EQ.  3))  THEN 
N=15 
MAX* 15 
NLESSK=5 

IF  (CGDE.EQ. 3)  THEN 
MAX* 15 
NLESSK=3 
ENDIF 

WRITE( 5 ,*) '  ' 

2  WRITE( 5 ,*) ' 2.  SELECT  BAUD  TYPE  (1-5)' 

READ( 5,*) BDTYPE 

IF  ((BDTYPE.  GT.  5).  OR.  (BDTYPE.  LT. 1))  THEN 
GO  TO  2 
ENDIF 
ENDIF 

IF  ( CODE. EQ. 6)  THEN 
MAX*  1 
ENDIF 


45 


o  o  o  nnn 


IF  ((CODE.  EQ.  2).  OR.  (CODE.  EQ.  4).  OR.  (CODE.  EQ.  5))  THEN 
N=31 
MAX=3 1 
KLESSK=5 

IF  (CODE. EQ. 4)  THEN 
MAX=31 
NLESSK=3 
END  IF 

WRITE( 5 ,*) '  ' 

221  WRITE (5 ,*) ' 2.  SELECT  BAUD  TYPE  (2-5)' 

READ (5 ,*)BDTYPE 

IF  ((BDTYPE.GT.  5).  OR.  (BDTYPE.LT.  2))  THEN 
GO  TO  221 
END  IF 
END  IF 

IF  ((CODE. GT.  0).  AND. (CODE. LT. 5))  THEN 
WRITE( 5 ,*) '  ' 

WRITE( 5 ,*) ' 3.  ERASE  CODEWORDS  WITH  UNCORRECTED  ERRORS  ?' 
WRITE ( 5  ,*) '  ENTER  "0"  TO  LEAVE  CODEWORDS  AS  IS' 

WRITE (5,*)’  ENTER  "l"  TO  ERASE’ 

RE AD (5,*) RUB 
END  IF 


WRITE( 5 ,*) 1  ' 

WRITE( 5 ,*) ' 4.  SELECT  SNR(DB,REAL  NUMBER)’ 
READ(5,*)SNR 


READ  IN  DATA  AND  CONSTANTS. . . 

CALL  BAUDS ( CODE , BDTYPE , K1 , K2 , PTS , CPB , M ) 

IF  ( CODE. EQ. 6)  THEN 
DO  108  1=0,3 

READ( 34 ,*) ( QPDIFF( I , J) , J=0 , 3 ) 

WR ITE ( 5 , * H  QPD I FF( I , J ) , J=0 , 3 ) 

READ( 14,*)QPRE( I) ,QPIM( I) 

QAMCOM( I )=CMPLX( QPRE( I ) ,QPIM( I ) ) 

108  CONTINUE 
GO  TO  597 
END  IF 

CALL  DATA( CODE ,N,GFA , TRACE , T4 , QAMCOM , QAMRE , QAMIM , QAMMAG , GFD1 6 , 
+  GFD32) 

597  NFIX=10*LOG10( (PTS/2)/(K2-Kl)/l. 0) 

DO  544  SNR=15. 0,22. 0,1. 0 

BYTERR=0 

N0ERR=0 

UNOERR=0 

TUOERR=0 

LRIERR=0 

FIXED=0 

DATERR=0 
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CWBERR=0 
DBERR=0 
ITER=0 
ERASE=0 
MEAN=0.  0 


DO  543  BDS=1 , 100 


DO 

15 

1=0, MAX-1 

CW(I)=0 

D( I )=0 

15 

CONTINUE 

DO 

16 

1=0 , PTS- 1 

TX(I)=(0.  0,0.  0) 

RX(I)=(0.  0,0.  0) 

T(I)=(0.  0,0.  0) 

16 

CONTINUE 

C  MAKE  UP  RANDOM  DATA,  ENCODE  IT,  LOAD  IT  IN  TRANSMITTER... 
IF  ((CODE. EQ. 0). OR. (CODE. EQ. 5))  THEN 

CALL  NOCODEC CPB , DO ,T , K1 , PTS , QAMCOM , DIN , N) 

GO  TO  797 
END  IF 

IF  (CODE. EQ. 6)  THEN 

CALL  QPSK(D0,K1,K2, PTS, QAMCOM, DIN, T) 

GO  TO  797 
END  IF 
PTR=K1 

DO  776  J=1,CPB 

CALL  R.NUN( MAX-NLESSK , D1 ) 

DO  1  1=0, MAX-NLESSK- 1 

D 1 ( I )=INT( (N+1)*D1( I ) ) 

1  CONTINUE 

CALL  INTRANC D 1 , D , MAX , NLESSK) 

CALL  ENCODE ( N , CV , D , GFA , CODE , NLESSK , MAX) 

DO  17  I=PTR , PTR+MAX- 1 

T( I )=QAMCOM(  CW( I -PTR) ) 

T( PTS-I )=CONJG(T( I ) ) 

DD( J , I -PTR)=D( I -PTR) 

C  C  WV  ( J ,  I  -  PTR )  =CW  ( I  -  PTR ) 

17  CONTINUE 

PTR=PTR+MAX 
776  CONTINUE 


TRANSMIT  THE  CODEWORD. . . 
797  CALL  BITREV( PTS ,M ,T ,TX) 
CALL  I NVFFT ( PTS , M , TX ) 


ANY  NOISE  IN  COMMUNICATION  CHANNEL?  (R(I)  IS  RECEIVED  VECTOR)... 
SUMSQ=0.  J 
SUM=0. 0  ' 

DO  13  1=0 , PTS- 1 
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TX( I )=TX( I)/2+C0NJG(TX( I) )/2 
SUM=SUM+REAL(TX( I ) ) 

SUMSQ=SUMSQ+( REAL( TX( I ) )**2 ) 

13  CONTINUE 

NVAR=SUMSQ/ ( PTS - 1 ) / ( 10**( ( SNR-NFIX) / 10.  0 ) ) 

CALL  GAUSS ( MEAN , NVAR , PTS , AWGN , TEMP 1 , TEMP2 ) 
DO  617  I=0,PTS-1 

TX( I )=TX( I )+AWGN( I ) 

617  CONTINUE 


RECEIVE  THE  DATA. . . 

CALL  B ITREV( PTS , M , TX , RX ) 

CALL  FFT(PTS ,M,RX) 

DECODE  EACH  RECEIVED  CODEWORD  SEQUENTIALLY.  .  . 

IF  (CODE. EQ. 0)  THEN 

CALL  WHAT16(RX, CPB , PTS ,K1 ,DOUT) 

PTR=K 1 
GO  TO  667 
END  IF 

IF  (CODE. EQ. 5)  THEN 

CALL  what 3 2 ( RX , N , PTS , QAMRE , QAMIM , DOUT , K 1 , CPB ) 
PTR=K1 
C-0  TO  668 
END  IF 

IF  (CODE. EQ. 6)  THEN 

CALL  VHATQP(RX,K] ,K2, PTS, DOUT) 

GO  TO  669 
END  IF 


PTR=K1 

DO  777  J=1 ,CPB 
B=0 
C=0 
K=0 

DO  881  1=0, MAX- 1 

CW(I)=CCWW(J,I) 

D( I )=DD( J , I ) 

DHAT( I )=0 
V(I)=0 
E( I )=0 
R( I )=0 

881  CONTINUE 

IF  ((CODE.  EQ.  1).  OR.  (CODE.  EQ.  3))  THEN 
CALL  WHAT 1 6 ( RX , MAX , PTS , PTR , R ) 

END  IF 

IF  ((CODE. EQ. 2). OR. (CODE.  EQ.  4))  THEN 

CALL  what 32 ( RX , MAX , PTS , QAMRE , QAMIM , R , PTR , MAX) 
END  IF 

PTR=PTR+MAX 


CALCULATE  THE  SYNDROME. . . 
ZEROS=0 
NEAkEN=0 
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00PS=0 

CALL  SYNDR0( R , S , N , NEAREN , OOPS , GFA , ZEROS , MAX , CODE ) 


CHECK  FOR  NO  ERRORS. . . 

199  IF  (ZEROS.  EQ.NLESSK)  THEN 

NOERR=NOERR+ 1  < 

CALL  RUB OUT ( RUB , R , V , DHAT , N , NLESSK , ZEROS , MAX) 
GO  TO  112 
END  IF 


GET  INTO  REST  OF  DECODING  ALGORITHM.  .  . 

IF  ((CODE.  EQ.  1).  OR.  (CODE.  EQ. 2))  THEN 
BAUDNO=J 
B=0 
C=0 
K=0 

CALL  DECTWO( R , S , V , DHAT , E .ERRORS , UNOERR , DUOERR , TRIERR 
+  , FIXED , B , C , K , ERASE , TRACE , N , GFA , GAMMA 1 , GAMMA 2 , GAMMA 3 

+  , RUB , NLESSK , ZEROS , CODE , MAX ) 

END  IF 

IF  ((CODE.  EQ.  3).  OR.  (CODE.  EQ.  4))  THEN 

CALL  DECONE (S.V, ZEROS , RUB , R , DHAT , N , NLESSK , MAX , GFA , CODE , 
+  UNOERR, DUOERR)  „ , 

END  IF 


NOW  CHECK  IF  ANY  SYMBOL  ERRORS  WERE  MADE  (COMPARE  ESTIMATE 
OF  CODEWORD  SENT  TO  ACTUAL  CODEWORD  SENT)... 

112  ERR0RS=0 

DO  61  1=0, MAX- 1 

IF  (V(I).NE.CW(I))  THEN 
B  YTE  RR=B  YTERR+ 1 
ERRORS=l 

IF  ((V(I).EQ.  -1).  AND.  ( (CODE.  EQ.  l).OR.  (CODE.  EQ.  3))) 

+  THEN 

CWBERR=MAX*4 
DBERR=(MAX-NLESSK)*4 
GO  TO  388 

END  IF 

IF  ((V(I). EQ.  -1).  AND.  ((CODE.  EQ.  2).  OR.  (CODE.  EQ.  4))) 

+  THEN 

CWBERR=MAX*5 
DBERR=(MAX-NLESSK)*5 
GO  TO  388 

ENDIF 

IF  ((CODE. EQ. 0). OR. (CODE. EQ. 1). OR. (CODE. EQ. 3))  THEN 
CWBERR=CWBERR+GFD16( V( I )  CW ( I ) ) 

GO  TO  133 
ENDIF 

CWBERR=CWBERR+GFD32(V( I ) ,CW( I ) ) 

133  ENDIF 

ERR0PS=1 

61  CONTINUE 

DO  336  1=0 ,MAX-NLESSK-1 
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231 

336 

C--- 

388 

C 

C113 


335 

99 

777 

667 


942 

668 


944 

669 


945 

543 


IF  (D(I).NE.DHAT(I))  THEN 
DATERR=DATERR+1 
IF  (DHAT(I). EQ. -1)  THEN 
DBERR=DBERR+4 
ELSE 

IF  ((CODE.  EQ.  0).  OR. (CODE. EQ. 1). OR. (CODE. EQ. 3))  THEN 
DBERR=DBERR+GFD 16 (D( I) , DHAT( I ) ) 

GO  TO  231 

END  IF 

DBERR=DBERR+GFD32( D( I ) ,DHAT( I ) ) 

END  IF 
END  IF 
CONTINUE 


BAUDNO=J 

IF  (ERRORS. EQ. 1)  THEN 

CALL  OUTPUT ( S , GAMMA 1 , GAMMA 2 , GAMMA3 .R.V.B.C.K.Ct^D, 
+  BAUDNO , DHAT , BDTYPE , MAX) 

END  IF 
NOP=0 

ITER=ITER+1 

CONTINUE 

IF  (CODE. EQ. 0)  THEN 
DO  942  1=0 , CPB- 1 

IF  (DIN(I).NE.  DOUT(I))  THEN 
DATERR=D ATERR+ 1 

DBERR=DBERR+GFD16(DIN( I) ,DOUT( I) ) 

END  IF 
CONTINUE 
END  IF 

IF  (CODE. EQ. 5)  THEN 
DO  944  1=0 ,CPB-1 

IF  (DIN(I).NE.  DOUT(I))  THEN 
DATERR=DATERR+ 1 

DBERR=DBERR+GFD32(DIN( I ) ,DOUT( I ) ) 

END  IF 
CONTINUE 
ENDIF 

IF  (CODE. EQ. 6)  THEN 
DO  945  1=0 ,K2-K1 

IF  (DIN(I).NE.  DOUT(I))  THEN 
DATERR=DATERR+ 1 

DBERR=DBERR+QPDIFF(DIN( I ) , DOUT( I ) ) 

ENDIF 

CONTINUE 

ENDIF 

CONTINUE 

IF  (CODE. EQ. 0)  THEN 

BE=CPB*4*(BDS-1)/1. 0 
BE=DBERR/BE 
WR I TE ( 3 3 , * ) SNR , BE 
ENDIF 

IF  (CODE. EQ. 5)  THEN 
BE=CPB''-'5*(  BDS - 1 ) / 1 .  0 
BE=DBERR/BE 
WRITE (33,*) SNR , BE 
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END  IF 

IF  (CODE. EQ. 6)  THEN 

BE=(K2-K1+1)*2*(BDS-1)/1.  0 
BE=DBERR/BE 
WR I TE ( 3  3 , * )  SNR ,  BE 
END  IF 

IF  ((CODE. EQ. 1). OR. (CODE. EQ. 3))  THEN 
BE=CPB*N*4*( BDS-1) / 1 .  0 
BE=DBERR/BE 
WRITE (33,*) SNR, BE 
END  IF 

IF  ( (CODE. EQ. 2). OR. (CODE. EQ. 4))  THEN 
BE=CPB*N*5*( BDS-1)/ 1.  0 
BE=DBERR/BE 
WRITE ( 33 ,*)SNR,BE 
END  IF 
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C756 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


CONTINUE 

WRITE( 6 ,*) '  SNR  (DB)  =  '  SNR 
WRITE( 6 ,*) '  BAUD  TYPE  2  ,BDTYPE 

IF  (CODE.  GT. 0)  THEN 
IF  (CODE. EQ. 1)  THEN 

WRITE(6,*) 'CODING  CHOICE  =  REED-SOLOMON  (15,10)' 
GO  TO  756 
END  IF 

IF  (CODE. EQ. 2)  THEN 

WRITE( 6,*) 'CODING  CHOICE  =  REED-SOLOMON  (31,26)' 
GO  TO  756 


END  IF 

IF  (CODE. EQ. 3)  THEN 

WRITE(6,*)' CODING  CHOICE  =  REED-SOLOMON  (15,12)' 

GO  TO  756 
END  IF 

IF  (CODE. EQ. 4)  THEN 

WKiTE(6,*) 'CODING  CHOICE  =  REED-SOLOMON  (31,28)' 

GO  TO  756 
END  IF 

WRITE(  6  ,*) '  CW  SYMBOL  ERRORS  AFTER  CORRECTION21'  ,BYTERR 
WRITE( 6 ,*) '  TOTAL  SYMBOLS  TRANSMITTED=' ,N*ITER 

WRITE ( 6 ,*) '  RESULTING  CODEWORD  BIT  ERRORS2' ,CWBERR 

IF  ((CODE. EQ. 1). OR. (CODE. EQ. 3))  THEN 

WRITE( 6 ,*) '  CODEWORD  BITS  TRANSMITTED2' ,N*ITER*4 

B  E=N* I TER*4 / 1 .  0 
BE=CWBERR/BE 
ELSE 

WRITE( 6 ,*) '  CODEWORD  BITS  TRANSMITTED=' ,N*ITER*5 

BE=N* I TER*5 / 1 .  0 
BE=CWBERR/BE 
END  IF 


WRITE( 6 ,*) '  CODEWORDS  WITH  NO  ERRORS=' ,NOERR 

WRITE(6 ,*) '  DETECTED  SINGLE  ERRORS=' .UNOERR 

WRITE ( 6 ,*) '  DETECTED  DOUBLE  ERRORS2' ,DUOERR 

IF  ((CODE. EQ. 1). OR. (CODE. EQ. 2))  THEN 

WRITE! 6 ,*) ’  DETECTED  TRIPLE  (OR  MORE)  ERRORS2' ,TRIERR 
WRITE ( 6 ,*) '  CORRECTED  TRIPLE  (BURST)  ERRORS2’ , FIXED 
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**v 


ENDIF 

WRITE( 6  ,*) '  ERASURES=' .ERASE 

WRITE (6,*) '  TOTAL  CODEWORDS  TRANSMITTED=’ , ITER 


ENDIF 

IF  (CODE. EQ. 0)  THEN 
WRITE(6,*)’ CODING 
VRITE( 6 ,*) ' 

WRITE( 6 ,*) ' 

WRITE( 6 ,*) ' 

WRITE( 6  ,*) ' 

ENDIF 

IF  (CODE. EQ. 5)  THEN 
WRITE( 6 ,*)' CODING 
WRITE( 6  ,*) ' 

WRITE( 6 ,*) * 

WRITE( 6 ,*) 1 
WRITE( 6  ,*) ' 

ENDIF 

WRITE( 6  ,*) '  ' 

WRITE (6,*)’ BIT  ERROR 
WRITE (6,*) ' 


CHOICE  =  UNCODED  (QAM-16)’ 

DATA  SYMBOL  ERRORS=' ,DATERR 
DATA  SYMBOLS  TRANSMITTED=' ,CPB*(BDS-1) 
DATA  BIT  ERR0RS=' ,DBERR 
DATA  BITS  TRANSMITTED= 1 , CPB*4*( BDS- 1 ) 


CHOICE  =  UNCODED  (QAM-32)’ 

DATA  SYMBOL  ERR0RS=' .DATERR 
DATA  SYMBOLS  TRANSMITTED=' ,CPB*(BDS-1) 
DATA  BIT  ERR0RS=' ,DBERR 
DATA  BITS  TRANSMITTED35'  ,CPB*5*(BDS-1) 


RATE  =  1 ,BE 

TOTAL  BAUDS  EXAMINED=’ ,BDS- 1 


r-.V 


WRITE( 33,*) ’TOTAL  BITS  IN  EACH  SNR=’ ,(K2-K1)*4*(BDS-1) 
ENDFILE( 6) 

CL0SE( 6) 

ENDFILE( 33) 

CLOSE( 33) 

END 


■  -.'r  iV  -.V  *;V  y?  y?  Vr  Vr  tV  Vr  Vr  V?  *  Vr  Vr  Vr  Vr  V?  -,V -V  Vr -.V  Vc  Vc  *V  Vc  *V  ‘V  Vr  Vr  Vc  Vc  Vr  Vc  Vc  Vc  Vc  Vr  Vc  Vr  Vr  Vc  -jV  “V  Vr  ■jV  Vr  Vr  Vc  Vc  Vr  Vr  ’jV  Vc  V  r  Vr  Vc  Vr  Vr 


SUBROUTINE  BAUDS(CODE , BDTYPE ,K1 ,K2 ,PTS ,CPB ,M) 

SUBROUTINE  TO  SET  SOME  INITIAL  VALUES,  DEPENDING  ON  BAUD  TYPE... 

CONSTANTS:  PTS . FFT  SIZE 

M . POWER  OF  2  (OF  PTS) 

K1 . FIRST  TONE  NUMBER  FOR  16-20KHZ  BAND  WITH 

SAMPLING  FREQ=61440HZ  AND  FFT  SIZE  PTS 

K2 . LAST  TONE  FOR  THIS  BAND 

CPB . CODEWORDS  PER  BAUD  (ONE  (15,10)  REED-SOLOMON 

CODEWORD  PER  15  TONES  IN  THIS  BAND) 

INTEGER  BDTYPE , K 1 , K2 , PTS , CPB , M , CODE 

IF  (BDTYPE.  EQ.  1)  THEN 
PTS=256 
M=8 
CPB=1 
Kl=68 
K2=83 
ENDIF 

IF  (BDTYPE. EQ. 2)  THEN 
PTS=5 12 
M=9 
CPB=2 
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Kl=135 

K2=166 

EN’DIF 

IF  ( BDTYPE. EQ.  3)  THEN 
PTS=1024 
M=10 
CPB=4 
Kl=269 
K2-332 
END  IF 

IF  (BDTYPE. EQ. 4)  THEN 
PTS=2048 
M-ll 
CPB=8 
K 1=5 3 7 
K2=664 
END  IF 

IF  (BDTYPE. EQ. 5)  THEN 
PTS=4096 
M=12 

IF  (CODE. EQ.  1)  THEN 
CPB=17 
ELSE 
CPB=16 
END  IF 
K1=1073 
K2=1328 
END  IF 

IF  ((CODE. EQ. 0). OR. (CODE. GE. 5))  THEN 
CPB=K2-K1+1 
END  IF 

IF  ((CODE. EQ. 2). OR.  (CODE.  EQ.  4))  THEN 
CPB=INT( CPB/2) 

END  IF 

RETURN 

END 


SUBROUTINE  PET:/  "OEE  ,N,GFA ,'rRACF.  ,74 , QAMCOM , QAMRE , QAMIM , OAMMAG , 
+  GFD16 , GFD32 ) 


THIS  SUBROUTINE  PROVIDES  DATA  TO  THE  PROGRAM  DEPENDING  ON  THE  CHOICE 
OF  CODE. 


ROWS . NUMBER  OF  ROWS  IN  DATA  FILE  FOR  EACH  ROW  IN  GF 

ADDITION  TABLE  (16  ENTRIES  IN  EACH  ROW  OF  DATA  FILE) 

QAMRE . REAL  PART  OF  A  QAM  CONSTELLATION  ELEMENT 

QAMIM .  IMAGINARY  PART  OF  A  QAM  CONSTELLATION  ELEMENT 

QAMCOM.  ...  THE  COMPLEX  ARRAY  OF  ( QAMRE , QAMRI ) 

QAMMAG. . . . MAGNITUDE  OF  THAT  QAM  ELEMENT 


INTEGER  N, SUM, ROWS, CODE 

INTEGER  GFA( 0: N,0: N) ,GFD16(0:  15,0:  15 ) ,TRACE( 0:  N) ,T4( 0:  N) 
INTEGER  GFD32(0: 31,0: 31) 

REAL  QAMRE ( 0: N) ,QAMIM(0: N) ,X,Y,QAMMAG(0: N) 
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COMPLEX  QAMCOM( 0:  N) 


IF  (CODE.  EQ.  1)  THEN 

READ( 7 ,*) (TRACE( I) , 1=0 ,N) 

READ( 7 ,*) (T4( I) , 1=0 ,N) 

ENDIF 

IF  (CODE. EQ. 2)  THEN 

READ( 12 ,*) (TRACE( I) , 1=0 ,N) 

ENDIF 

ROWS=INT(N+l)/16 
DO  10  1=0, N 
SUM=0 

DO  20  J=1 .ROWS 

IF  ((CODE.  EQ.  0).OR.  (CODE.  EQ.  1).  OR.  (CODE.EQ.  3))  THEN 
READ( 9 ,*) (GFA( I ,K) ,K=0,15) 

READ( 1 ,*) (GFD16( I ,K) ,K=0 , 15) 

GO  TO  20 
ENDIF 

RE AD ( 1 1 , * )  ( GFA  ( I , K ) , K=SUM , SUM+ 1 5 ) 

READ( 31,*)  ( GFD32 ( I , K) , K=SUM , SUM+15 ) 

SUM=SUM+16 
20  CONTINUE 

10  CONTINUE 

DO  40  1=0, N 

IF  ((CODE.EQ.  0).  OR.  (CODE.EQ.  1).  OR.  (CODE.EQ.  3))  THEN 
READ( 8 , * ) QAMRE (I) ,QAMIM(I) , QAMMAG( I ) 

GO  TO  35 
ENDIF 

READ( 13,*) QAMRE ( I ) , QAM I M ( I ) , Q AMMAG( I ) 

35  X=QAMRE( I) 

Y=QAMIM( I) 

QAMCOM( I )=CMPLX(X , Y) 

40  CONTINUE 
RETURN 
END 


0*Vf****-.V**Vf****Vf***Vf**Vf******'!V**Vf*,>V*****?V')V'sV'>V***iV*,>V**?V**'iV,»V')V****Vc*iViWfiV'(Wf»V 


SUBROUTINE  NOCODE( CPB , DO ,T,K1 , PTS , QAMCOM , DIN , N) 


THIS  SUBROUTINE  GIVES  RANDOM  DATA  FOR  THE  UNCODED  CASES  AND 
LOADS  IT  INTO  THE  TRANSMITTER.  .  . 


DIN . INTEGER  REPRESENTING  AN  ELEMENT  OF  GF(2**M) 

DO . REAL  NUMBER  EQUAL  TO  INTEGER  DIN  (FOR  INV-FFT) 

T . CONJUGATE  SYMMETRIC  SPECTRUM  FOR  INV-FFT 


INTEGER  CPB ,K1 ,DIN( 0:  255) , PTS ,N 
REAL  D0( 0: 255 ) 

COMPLEX  T(0: 4095) ,QAMCOM( 0: 31) 

CALL  RNUN(CPB.DO) 

DO  10  1=0 ,CPB-1 

D0( I j=INT( (N+1)*D0( I) ) 

DO  20  J=0,N 

IF  (D0(I).  EQ.  J/l.  0)  THEN 
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DIN(I)=J 
END  IF 

20  CONTINUE 

T(K1+I)=QAMC0M(DIN( I) ) 

T( PTS-K1 - I )=C0NJG(T( Kl+I ) ) 
10  CONTINUE 
RETURN 
END 


Q  Vc V?  “.'f  V?  V:  V?  V?  V?  “V  V?  *  Vr  Vr  Vr  Vr  V?  Vr  Vr  V?  V:  “V  Vr  V?  V?  *  i:  Vr  '.V  V:  ‘V  -V  *  Vf  V?  Y?  Yc  Vr  V?  Yr  Yr  Yr  Yr  *  Vr  Yc  Yr  Yr  Yr  Yr  Yc  Yr  Yr  Yr  Yc  Yc  Yr  Yc  Yr  *  Yr  *  Yc  Yr  Yr  Yr  *  Yc  Yc  * 

SUBROUTINE  QPSK( DO ,K1 ,K2 , PTS , QAMCOM , DIN , T) 


THIS  SUBROUTINE  GIVES  RANDOM  DATA  FOR  UNCODED  QPSK  AND 
LOADS  IT  INTO  THE  TRANSMITTER.  .  . 


DIN . INTEGER  REPRESENTING  AN  ELEMENT  OF  GF(2**M) 

DO . REAL  NUMBER  EQUAL  TO  INTEGER  DIN  (FOR  INV-FFT) 

T . CONJUGATE  SYMMETRIC  SPECTRUM  FOR  INV-FFT 


INTEGER  K1,K2,DIN(0: 255) ,PTS 
REAL  DO ( 0 : 255 ) 

COMPLEX  T( 0 : 4095) ,QAMC0M(0:  31) 

CALL  RNUN( K2-K1+1 .DO) 

DO  10  1=0 ,K2-K1 

D0( I )=INT(4*D0( I) ) 

DO  20  J=0 , 3 

IF  (D0( I).  EQ.  J/l.  0)  THEN 
DIN( I)=J 
ENDIF 

20  CONTINUE 

T(K1+I )=QAMCOM(DIN(  I ) ) 

T( PTS-K1-I )=CONJG(T( Kl+I ) ) 
10  CONTINUE 
RETURN- 
END 


*c*  r1;  r  *.  CYC"  ric“  a 


■  Yr  Yr  Yc  Yr  -V  Yc  Yr  Yr  Yr  Yr  Vr  Yr  Yc  Yc  Yc  '/V  Yr  “V  “V  Yc  Yr 


YcYr  Yc  YcYr  *  Yr  *  Yc  “V  Yc  Yr  *  YrYrYr  Yr  Yr  Yr -V  YcYr  Yr  * 


SUBROUTINE  ENCODE (N ,  CW,D,GFA, CODE , NLESSK ,MAX) 

THIS  SUBROUTINE  IS  THE  REED-SOLOMON  ENCODER  (SYSTEMATIC  FORM). 


G0-G4 . COEFFICIENTS  OF  GENERATOR  POLYNOMIAL 

SR . SHIFT  REGISTERS 


INTEGER  N , INPUT ,00,01,02,03,04, CODE , NLESSK , MAX 
INTEGER  CW( 0 : MAX - 1 ) ,D(0:  MAX-1) ,GFA(0:N,0:  N) ,SR(5) 

C  SET  GENERATOR  POLYNOMIAL  CONSTANTS... 

Gl=5 


c 


G2=12 

G3=12 


G4=5 

IF  ( CODE. EQ. 2)  THEN 
Gl=14 
G2=18 
G3=18 
G4=14 

END  IF 

IF  (CODE. EQ. 3)  THEN 
G0=4 
Gl=12 
G2=l  1 

END  IF 

IF  (CODE. EQ. 4)  THEN 
G0=4 
Gl=13 
G2=12 


i 


END  IF 

ENCODE  THE  DATA  FOR  THE  SYSTEMATIC  (N,N-5)  CODE... 


411  DO  8  J=1 .NLESSK 
SR( J)=0 

8  CONTINUE 

DO  9  J=0 ,MAX-NLESSK- 1 

CW ( MAX - 1 - J ) =D ( MAX -NLESSK - 1  - J ) 

INPUT=GFA( SR( 1 ) ,D(MAX-NLESSK-1-J) ) 

IF  (INPUT.  EQ.  O')  THEN 
DO  40  K=l, NLESSK- 1 
SR( K)=SR( K+l ) 

40  CONTINUE 

SR(NLESSK)=I NPUT 
GO  TO  9 
END  IF 

IF  (NLESSK. EQ. 5)  THEN 

SR( 1)=GFA( SR( 2) ,MOD( INPUT+G4-2 ,N)+1) 
SR( 2)=GFA( SR( 3) ,MOD( INPUT+G3-2 ,N)+1) 
SR(3)=GFA(SR(4) ,MOD( INPUT+G2-2 ,N)+1) 
SR( 4)=GFA( SR( 5 ) ,MOD( INPUT+G1-2 ,N)+1) 
SR( 5 )=INPUT 
END  IF 

IF  (NLESSK. EQ. 3)  THEN 

SR( 1)=GFA( SR( 2) ,MOD( INPUT+G2-2 ,N)+1) 
SR( 2 )=GFA( SR( 3 ) , MOD( INPUT+G1-2 ,N)+1 ) 
SR( 3)=M0D( INPUT+GO-2 ,N)+1 
END  IF 

9  CONTINUE 

DO  11  1=0 , NLESSK- 1 
CW( I )=SR( NLESSK- I ) 

11  CONTINUE 

RETURN 

END 

'  w  .  fm  ft*#  »*•  k*.  J.  J.  w  J.  k'k  k'k  y.  k'k  J.  J-  J-  kU  Jk  j-  Jk  k'k  kU  J* 

SUBROUTINE  WHAT16(X,MAX,L,PTR,R) 
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SUBROUTINE  TO  DECIDE  WHICH  ELEMENT  OF  GF(N+1)  EACH  ELEMENT 
OF  THE  RECEIVED  DATA  CORRESPONDS  TO. 

PTR . CURRENT  CODEWORD  STARTING  TONE 

X() . AN  ELEMENT  FROM  THE  RECEIVER  (FFT) 

MAX . CODEWORD  BLOCK  LENGTH 

L . FFT  SIZE 

R() . THE  RECEIVED  CODEWORD  IN  TERMS  OF  ELEMENTS  OF  GF(N+1) 


INTEGER  MAX,R(0:  MAX-1) ,PTR 
COMPLEX  X( 0: L-l) 

REAL  A , B , Z ,MID 

MID=150. 0 

DO  10  J=PTR , PTR+MAX- 1 
A=REAL( X( J) ) 

B=AIMAG(X( J) ) 

Z=CABS(X( J) ) 

IF  ((A. GE. 0). AND. (B. GE. 0))  THEN 

IF  ( ( Z. GE. MID). AND. (A. GE. B))  THEN 
R( J-PTR)=4 
GO  TO  10 
END  IF 

IF  ( ( Z. GE. MID). AND.  (A.  LT.  B) )  THEN 
R(J-PTR)=7 
GO  TO  10 
END  IF 

IF  ((Z. LT. MID).  AND. (A.  GE. B))  THEN 
R( J-PTR)=0 
GO  TO  10 
END  IF 

R( J-PTR)=3 
GO  TO  10 
END  IF 

IF  ((A.  GE.  0).  AND.  (B.  LT.  0))  THEN 

IF  ( ( Z . GE. MID). AND. (A. GE. ABS(B)))  THEN 
R( J-PTR)=15 
GO  TO  10 
END  IF 

IF  ( (Z. GE. MID). AND. (A.  LT. ABS(B) ) )  THEN 
R(J-PTR)=14 
GO  TO  10 
END  IF 

IF  ((Z.  LT.  MID).  AND.  (A. GE.  ABS(B)))  THEN 
R( J-PTR)=1 
GO  TO  10 
ENDIF 

R( J-PTR)=9 
GO  TO  10 
ENDIF 

IF  ((A. LT. 0).  AND.  (B.  LT.  0))  THEN 

IF  ((Z.GE.MID). AND. (ABS(A). GE. AES(B) ))  THEN 
R( J-PTR)=8 
GO  TO  10 
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ENDIF 

IF  ((Z.GE.MID).AND.  (ABS(A).LT.  ABS(B)))  THEN 
R( J-PTR)=13 
GO  TO  10 
ENDIF 

IF  (CZ.LT.MID).  AND.  (ABS(A).GE.  ABS(B)))  THEN 
R(J-PTR)=5 
GO  TO  10 
ENDIF 

R( J-PTR)=1 1 
GO  TO  10 
ENDIF 

IF  ((A.  LT.  0).  AND.  (B.  GE.  0))  THEN 

IF  ( (  Z.  GE.  MID) .  AND.  (ABS(A).GE.B))  THEN 
R(J-PTR)=10 
GO  TO  10 
ENDIF 

IF  ((Z.  GE.  MID).  AND.  (ABS(A). LT.  B))  THEN 
R( J-PTR)=12 
GO  TO  10 
ENDIF 

IF  ( ( Z.  LT.  MID).  AND.  (ABS(A). GE. B) )  THEN 
R( J-PTR)=2 
GO  TO  10 
ENDIF 

R( J-PTR)=6 
ENDIF 
CONTINUE 

RETURN 

END 


*  V  Vr  *V  * V  Vr  Vr  -A-  -,Y  Vc  Vc  ‘.V  Vc  Vc  Vc  -V  Vc  Vc  Vc  Vc  Vc  Vc  V r  Vc  Vc  Vc  Vc  Vc  Vc  Vc  Vc  Vc  *  Vc  ‘V  Vc  Vc  Vc  Vc  Vc  Vc  Vc  Vc  Vc  *  Vc  ‘.V  Vc  Vc  Vc  Vc  V:  Vc  V;  Vc  Vc  Vc  Vc  Vc  Vc  Vc  Vc  Vc  Vc  ‘V  V c  Vc 


SUBROUTINE  WHATQP( X ,K1 ,K2 , L,R) 


SUBROUTINE  TO  DECIDE  WHICH  ELEMENT  OF  GF(N+1)  EACH  ELEMENT 
OF  THE  RECEIVED  DATA  CORRESPONDS  TO. 


X() . AN  ELEMENT  FROM  THE  RECEIVER  (FFT) 

K1 . FIRST  TONE  IN  BAUD 

K2 . SECOND  TONE  IN  BAUD 

L . FFT  SIZE 

R( ) . THE  RECEIVED  CODEWORD  IN  TERMS  OF  ELEMENTS  OF  GF(N+1) 


INTEGER  R(0: K2-K1+1) 

COMPLEX  X( 0: L- 1) 

REAL  A , B 

DO  10  J=K1 ,K2 

A=REAL(X( J) ) 

B=AIMAG(X( J) ) 

IF  ((A. GE. 0). AND. (B.  GE. 0))  THEN 
R(J-K1)=0 
GO  TO  10 
ENDIF 

IF  ( f  A.  LE.  0). AND. ( B. GE. 0) )  THEN 
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R(  J-K1)  =  1 
GO  TO  10 
END  IF 

IF  ((A.  LE.  0).  AND.  (B.  LE.  0))  THEN 
R( J-K1)=2 
GO  TO  10 
END  IF 
R( J-Kl )=3 
10  CONTINUE 
RETURN 
END 


Q  **VrsV*^V* 


Vr  ■>’?  -V  V?  Vc  Vf  V?  “V  *'r  *  ■ 


Vr*  VrVrVrVr  iWrVr-.V  ^WfVrs’oWrVoV  VoWc  *  ***■>>**  yrVc*  * 


SUBROUTINE  VHAT32 ( X ,  N ,  L , QAMRE , QAMIM , R , PTR , MAX) 


SUBROUTINE  TO  FIND  CLOSEST  QAM  CONSTELLATION  ELEMENT  BY  SUBTRACTING 
THE  MAGNITUDE  OF  THE  RECEIVED  COMPLEX  SIGNAL  FROM  THE  PRODUCT  OF 
THE  COMPLEX  CONJUGATE  OF  THE  RECEIVED  SIGNAL  AND  EACH  ELEMENT  OF 
THE  QAM  CONSTELLATION,  AND  TAKING  THAT  CONSTELLATION  ELEMENT 
WITH  THE  SMALLEST  RESULT  AS  AN  ESTIMATE  OF  THE  RECEIVED  SIGNAL. 
(ONLY  USED  WITH  32 -QAM) 


INTEGER  N , L , K( 0:  MAX- 1 ) , WHO , PTR 
COMPLEX  X( 0:  L-l) 

REAL  NEW , A , NEWR , NEW I , QAMRE ( 0 :  N) ,QAMIM(0:  N) 

DO  10  I =PTR , PTR+MAX - 1 

NEWR=ABS( QAMREC  0) -REAL(X( I ) ) )** 2 
NEWI=ABS(QAMIM( 0) -AIMAG(X( I) ) )**2 
NEW=SQRT( NEWR+NEWI ) 

WH0=0 

DO  20  J=1 ,N 

A=ABS( QAMREC  J ) -REAL( X( I ) ) )**2 
A=A  +  ABS( QAMIM( J) -AIMAG(X( I) ) )**2 
IF  (SQRT(A).  LT.  NEW)  THEN 
NEW=  A 
WHO=J 
END  IF 

20  CONTINUE 

R(  I -PTR)=WKO 
10  CONTINUE 
RETURN 
END 


»V  *  :V?V  ?V  *  ft : 


Vc'VrVrVr^V^V^V 


SUBROUTINE  SYNDROC R , S , N , NEAREN , OOPS , GFA , ZEROS , MAX , CODE ) 


THIS  SUBROUTINE  COMPUTES  THE  SYNDROME  VECTOR  (FOR  DECODING) 

S . SYNDROME  VECTOR 

ZEROS . NUMT  _  1  OF  SYNDROME  ELEMENTS  EQUAL  TO  ZERO 

NEAREN . FLAG  FROM  DECODING  SUBROUTINE,  =1  IF  THIS  IS  THE 

SECOND  SYNDROME  COMPUTATION  FOR  A  GIVEN  CODEWORD 
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(RECHECKING  CORRECTIONS  HADE  FOR  HIGHER  ORDER  ERRORS) 

OOPS . FLAG  TO  DECODING  SYNDROME  THAT  AN  ELEMENT  OF  SYNDROME 

IS  NOT  EQUAL  TO  ZERO  (WHILE  PERFORMING  SECOND  SYNDROME 
COMPUTATION) 

LOW . FIRST  ELEMENT  OF  SYNDROME  (DIFFERS  BETWEEN  LENGTH  31 

AND  LENGTH  15  CODES) 


INTEGER  N , ZEROS , NEAREN , OOPS , X , ALPHA , TEMP , SUM , CODE 
INTEGER  R(0: MAX-1) ,S( -2:  2) ,GFA( 0:  N , 0:  N) , LOW 
OOPS=0 
DO  5  I=-2 , 2 
S( I )=0 

5  CONTINUE 
ZEROS=0 

IF  (CODE. GT. 2)  THEN 
LOW=0 
ELSE 
LOW=-2 
END  IF 

DO  20  ALPHA=LOW , 2 
SUM=0 

DO  30  1=0, MAX- 1 

IF  (R(I).EQ.  0)  THEN 
GO  TO  30 
END  IF 
TEMP=0 

IF  (ALPHA.  LT.  0)  THEN 
TEMP=ALPHA+N 
X=R( I ) - 1+TEMP* I 

try  cr 

X=R( I ) - 1+ALPHA" I 
END  IF 
X=MOD(X,N) 

X=X+1 

SUM=GFA( SUM ,X) 

30  CONTINUE 

S(ALrHA)=SUM 

IF  ((SUM. NE.O). AND.  (NEAREN.  EQ.  1))  THEN 
OOPS=l 
GO  TO  21 
END  IF 

IF  (S(ALPHA).  EQ.  0)  THEN 
ZEROS=ZEROS+l 
END  IF 

20  CONTINUE 

21  RETURN 
END 


:V  Vr  iV  *V  Vr  -V Vr  'V  Vr  V:  Vr  ‘V  Vr  Vr 


VcVrVr  VrVr^VVrVrVrVrVrVrVrVrVfVfVrVrVr^VyrVrVrVrVrVrVr  VrVr  Vr  Vr  Vr  Vr  Vr  Vr  V-VoY  ?Y  Vr  -jV  Vr  ?Y  tY  r’r  Vr  ?Y  Vr  Vr 


SUBROUTINE  DECTWO(R , S , V , DHAT,E , ERRORS ,UN0ERR ,DUOERR ,TRIERR 
+  .FIXED , B ,C ,K, ERASE , TRACE , N , GFA , GAMMA 1 , GAMMA 2 , GAMMA3 

+  , RUB, NLESSK, ZEROS, CODE, MAX) 


DECODING  SUBROUTINE  FOR  DEC/TED  CODES 
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R . RECEIVED  CODEWORD  VECTOR 

S . SYNDROME  VECTOR 

E . ESTIMATED  ERROR  VECTOR 

V . ESTIMATED  CODEWORD  VECTOR 

DHAT . ESTIMATED  DATA  VECTOR  (FROM  V) 

ERRORS . 17 LAG  FOR  PARAMETERS  OUTPUT  IN  CASE  THERE  ARE 

DECODING  ERRORS  (ERRORS=l) 

UNOERR . SINGLE  ERRORS 

DUOERR . DOUBLE  ERRORS 

TRIERR . UNCORRECTABLE  TRIPLE  ERRORS 

FIXED . CORRECTED  TRIPLE  BURST  ERRORS 

B,C,K . PARAMETERS  USED  IN  DOUBLE  ERROR  CORRECTION 

GAMMA1-3 . PARAMETERS  USED  IN  DOUBLE  ERROR  CORRECTION 

TRACE . TRACE  OF  AN  ELEMENT  OF  GF(2**M) 

ERASE . NUMBER  OF  ERASURES  (IF  ERASING  TRIPLE-ERRORS  IS 

CHOSEN  (IF  RUB=1) 


INTEGER  GAMMA 1 , GAMMA2 , GAMMA3 , UNOERR , DUOERR , TRIERR , N , MAX 
INTEGER  FIXED, NOP, B,C,K,V1,V2,W3, XI, X2 ,X3 ,Y1 , Y2 ,Y3 , TEST 
INTEGER  ROOT( 2) , TRACE! 0: N) ,GFA(0:  N,0:  N) 

INTEGER  R( 0:  N-l)  ,S(-2:  2),V(0:  N-l)  ,DHAI'(0:  N-l)  ,E(0:  N-l) 
INTEGER  Z3 , W , X , Y , ERASE , A 1 , A2 , ERRORS , OOPS , ZEROS , RUB , CODE 


C  CHECK  FOR  3  OR  MORE  ERRORS  FROM  SYNDROME.  .  . 

IF  (ZEROS. GT. 2)  THEN 
TRIERR=TRIERR+1 

CALL  RUB0UT( RUB , R , V , DHAT , NLESSK , ZEROS , MAX ) 
GO  TO  333 
END  IF 


36  NOP=0 

IF  (S( 1)*S( -2). EQ. 0)  THEN 
A1=0 
ELSE 

Al=MOD( S( 1)+S( -2) -2 ,N)+1 
END  IF 

IF  (S( -1)"S(0). EQ. 0)  THEN 
A2=0 
ELSE 

A2=MOD(S(-1)+S(0)-2,N)+1 
END  IF 

GAMMA 1=GFA(A1 ,A2) 

IF  ((S(2)*S(-2)).EQ. 0)  THEN 
A1=0 
ELSE 

A1=M0D( S( 2)+S(-2)-2,N)+l 
ENDIF 

IF  (S(0).  El.  0)  THEN 
A2 =0 
ELSE 
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A2=M0D(S(0)*2-2,N)+1 
END  IF 

GAMMA2=GFA(A1,A2) 

IF  ( ( S( 0)*S( 1) ) . EQ. 0)  THEN 
A1=0 
ELSE 

A1=M0D(S(0)+S(1)-2,N)+1 
END  IF 

IF  ( ( S( 2)*S( -1) ). EQ. 0)  THEN 
A2=0 
ELSE 

A2=M0D( S( 2)+S( -1)-2,N)+1 
ENDIF 

GAMMA3=GFA( A1 , A2 ) 


STEP# 4:  CHECK  FOR  A  SINGLE  ERROR.  .  . 


IF  (( GAMMA  1.  E0.  0).  AND.  (GAMMA2.  EQ.  0).  AND.  (GAMMA 3.  EQ.  0))  THEN 
Y=M0D(S( 1)+N-S(0) ,N) 

E(Y)=S(0) 

U.N0ERR=UN0ERR+1 
GO  TO  59 
ENDIF 


STEP??5:  CHECK  FOR  3  OR  MORE  ERRORS.  .  . 

IF  ( ( ( GAMMA 1 . NE. 0). OR.  (GAMMA2. NE.  0). OR. (GAMMA3. NE.  0)).  AND.  ( GAMMA 1 
+  *GAMMA2'VGAMMA3.  EQ.  0)  )  THEN 

TRIERK=TRIERR+1 
ERASE=ERASE+1 

CALL  RUE  OUT ( RUB  ,  R ,  V ,  DHAT ,  N'LESSK ,  ZEROS ,  MAX) 

GO  TO  333 
ENDIF 


STEP?- 6:  CALCULATE  CONSTANTS  IN  QUADRATIC  FORMULA.  .  . 

IF  (GAMMA2.EQ. 0)  THEN 
B=0 
ELSE 

B=M0D( GAMMA2+N-GAMMA1 ,N)+1 
ENDIF 

IF  (GAMMA3. EQ. 0)  THEN 
C=0 
ELSE 

C=M0D ( GAMMA3+N - GAMMA 1 , N ) + 1 
ENDIF 

IF  (C.  EQ.  0)  THEN 
K=0 

GO  TO  37 
ENDIF 

IF  (B. EQ. 1)  THEN 
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K=C 

GO  TO  37 
END  IF 

K=MOD( C - l+3*N-2*( B - 1 ) ,N)+1 
37  NOP=0 

IF  ( (TRACE(K).  EQ. 1).  AND.  (K. ME. 6) )  THEN 
TRIERR=TRIERR+1  f 

ERASE=ERASE+1 

CALL  RUB OUT ( RUB ,R , V , DHAT , NLESSK , ZEROS , MAX) 
GO  TO  333 
END  IF 


LOOK  FOR  A  CORRECTABLE  TRIPLE  BURST  ERROR.  .  . 


IF  (K.  EQ.  6)  THEN 
W1=S(-1)-1 
W2=S(1)-1 
W3=M0D(W1+W2,N) 

X1=S( -2) -1 
X2=S(2)-1 
X3=MOD(Xl+X2 ,N) 

Y1=S( 0) - 1 
IF  (N. EQ. 15)  THEN 
Y2=M0D(2*Y1+3,N) 
Y3=M0D(2*Y1+6,N) 

ELSE 

Y2=M0D(2*Y1+20,N) 

Y  3=MOD ( 2* Y 1+2  3 , N ) 

END  IF 

IF  ( (W3.  EQ.  Y2).  AND.  (X3.  EQ.  Y3) )  THEN 
F I XED=F I XE  D+ 1 
IF  (N. EQ. 15)  THEN 

Z  3=M0D ( V2 - Y 1+3*N - 1 0 , N ) 

ELSE 

Z  3=M0D( V2 - Y 1+3*N - 1 1 , N) 

END  IF 
Ef Z3)=S(0) 

E( Z3+1 )=S( 0) 

Ef  Z3+2)=S( 0) 

GO  TO  59 
END  IF 


END  IF 


FIND  THE  ERRORS  USING  THE  QUADRATIC  FORMULA.  .  . 

TEST=0 
ROOT( 1 )=0 
ROOT( 2)=0 
DO  50  1=1, N 

W=MOD( 2* ( I - 1 ) , N) + 1 
IF  (B. EQ. 1)  THEN 
X=I 

GO  TO  41 
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END  IF 

IF  (B.EQ.  0)  THEN 
X=0 

GO  TO  41 
END  IF 

X=M0D(B+I-2+N,N)+l 

41  Y=GFA(W,X) 

Z=GFA( Y,C) 

IF  ((Z.EQ.  0).  AND.  (TEST.  EQ.O))  THEN 
R00T( 1)=I - 1 

42  TEST=1 
GO  TO  50 

END  IF 

IF  ((Z.EQ.  0).  AND.  (TEST.  EQ.  1))  THEN 
R00T(2)=I-1 
43  GO  TO  55 

END  IF 
50  CONTINUE 


55  NOP=0 


C  ROOT(l)  AND  ROOT( 21  GIVE  DOUBLE  ERROR  LOCATIONS,  NOW  FIND  MAGNITUDE 
IF  (S(0). EQ- 0)  THEN 
Al=0 
ELSE 

A1=M0D( S( 0) -1+R00T( 2) ,N)+1 
END  IF 

IF  (S(l).EQ.  0)  THEN 
A2=A1 
ELSE 

A2=GFA( A1 ,S( 1) ) 

END  IF 

X=M0D( A2+N-B ,N)+1 

E(ROOT( 1) )=X 

E( ROOT( 2) )=GFA( S( 0) ,X) 

IF  (E(R00T(1))*E(R00T(2)).EQ.  0)  THEN 
TR I ERR=TR I ERR+ 1 
ERASE=EK.iSE+l 

CALL  RUBOUT( RUB , R , V , DHAT , NLESSK , ZEROS , MAX) 

GO  TO  333 
END  IF 

DUOERR=DUOERR+ 1 
59  NOP=0 


STEP# 7:  NOW  ADD  ESTIMATED  ERROR  VECTOR  TO  RECEIVED  VECTOR. . . 
ERRORS=0 

72  DO  60  1=0, N-l 

V( I)=GFA(R( I) ,E( I) ) 

IF  (I.GT.  4)  THEN 
DHAT( I -5 )=V( I ) 

END  IF 
60  CONTINUE 
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c . 

C  STEP,-!-1 8:  MAKE  LAST  CHECK  FOR  3  OR  MORE  ERRORS.  .  . 


NEAREN=1 

ZEROS=0 

00PS=G 

CALL  SYNDR0( V , S , N , NEAREN , OOPS , GFA , ZEROS , MAX , CODE ) 
IF  (OOPS.EQ. 1)  THEN 
TRIERR=TRIERR+1 
ERASE=ERASE+1 
DUOERR=DUOERR- 1 

CALL  RUBOUTC RUB , R , V , DHAT , NLESSK, ZEROS , MAX) 

GO  TO  333 
END  IF 


333 


C***- 


RETURN 

END 

■  Vr  Vr;V  Vr  Vc  ■: Wf  Vr  Vc  VrVoV  •) WrVr -sV  VfV?  Vr*  VcVr  ■jV  ^VV?  Vf  Vr  ■ 


■VoWr 


*•>'?  VoWr  Vc  Vc  Vr  Vr  *  *  Vc  iV  *Wc  *  Vr*  *WnV  ■sV  Ve  *  iV-;V  Vc  Vc  Vr 


SUBROUTINE  DECONE ( S , V , ZEROS , RUB , R , DHAT , N , NLESSK , MAX , GFA , CODE , 
+  UNOERR , DUCERR ) 


DECODING  SUBROUTINE  FOR  SEC/DED  CODES 


R . RECEIVED  CODEWORD  VECTOR 

S . SYNDROME  VECTOR 

V . ESTIMATED  CODEWORD  VECTOR 

DHAT . ESTIMATED  DATA  VECTOR  (FROM  V) 

DECODING  ERRORS  (ERR0RS=1) 
UNOERR . SINGLE  ERRORS 


DUOERR . DOUBLE  ERRORS 

CHOSEN  (IF  RUB=1) 


INTEGER  ZEROS , N , DUOERR ,X, Y, UNOERR , MAX , CODE , RUB , NLESSK 

INTEGER  S( -2  ) ,GFA(0:  N,0:  N) ,V(0:  MAX-1) ,DHAT( 0: MAX-1) ,R(0: MAX-1) 

INTEGER  NEAREN, OOPS 

IF  (ZEROS. EQ. 1)  THEN 
DUOERR=DUOERR+l 

CALL  RUB OUT ( RUB , R , V , DHAT , NLESSK , ZEROS , MAX ) 

GO  TO  150 
END  IF 
S( -2)=0 
S(-1)=0 

DO  10  1=0, MAX- 1 
V( I )=R( I ) 

10  CONTINUE 

IF  ((S(O).NE. 0). AND. ( ZEROS. EQ. 2))  THEN 
V( 0 )=GFA( S( 0) ,R(0) ) 

UNOERR=UNOERR+ 1 
GO  TO  100 
END  IF 

IF  ( ( S( 1) .  NE.  0) .  AND.  ( ZEROS. EQ.  2) )  THEN 
V( 1 )=GFA( S( 1 ) ,R( 1) ) 

UNOERR=UNOERR+ 1 


65 


o  o  o  o 


GO  TO  100 
END  IF 

IF  ( ( S ( 2 ) .  NE.  0). AND. (ZEROS. EQ. 2))  THEN 
V( 2)=GFA( S( 2) ,R(2) ) 

UNOERR=UNOERR+l 
GO  TO  100 
END  IF 

X=MOD(S( 1)+N-S(0) ,N) 

Y=MOD(S(2)+N-S( 1) ,N) 

IF  (X. EQ. Y)  THEN 
C  x=X+3 

V(X)=GFA(R(X) , S(0)) 

UN  OERR=UNOERR+ 1 
GO  TO  100 
ELSE 


DUOERR=DUOERR+l 

CALL  RUBOUT( RUB , R , V , DHAT , NLESSK , ZEROS , MAX) 

GO  TO  150 
END  IF 

C  LAST  CHANCE  TO  CHECK  FOR  2  OR  MORE  ERRORS.  .  . 

100  NEAREN=1 
ZEROS=0 
OOPS=0 

CALL  SYNDRO( V ,  S  ,  N , NEAREN , OOPS , GFA , ZEROS , MAX , CODE ) 
IF  (OOPS. EQ. 1)  THEN 
DUOERR=DUOERR+ 1 
UNOERR=l'NOERR- 1 

CALL  RUBOUTC RUB , R , V , DHAT , NLESSK , ZEROS , MAX) 

GO  TO  150 
END  IF 

DO  20  1=3, MAX- 1 

DHAT( I -3)=V( I ) 

20  CONTINUE 
150  RETURN 
END 


C 


Msir 


*,V  Vr  Vc  *  Mr  *  Mn't  MriT  Vc  ie  Vc  ic  Mr  V?  *  Me  “V 


SUBROUTINE  RUBOUT( RUB , R , V , DHAT , NLESSK , ZEROS , MAX) 


THIS  SUBROUTINE  ERASES  CODEWORDS  WITH  MORE  ERRORS  THAN  A  GIVEN  CODE 
CAN  CORRECT. 


INTEGER  NLESSK, ZEROS, MAX 

INTEGER  RUB ,R(0: MAX-1) ,V(0:  MAX-1) ,DHAT(0:  MAX-1) 

IF  ((RUB.  EQ.  0).  OR.  (ZEROS.  EQ.  NLESSK))  THEN 
DO  14  1=0, MAX- 1 
V( I)=R( I) 

IF  (I.  GT. NLESSK- 1)  THEN 
DHAT( I -NLESSK)=R( I ) 

ENDIF 

14  CONTINUE 

GO  TO  50 
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END  IF 

IF  (RUB.  EQ. 1)  THEN 
DO  101  1=0 , MAX-1 
V(I)=-1 

IF  (I.  LE.MAX-NLESSK-1)  THEN 
DHAT( I )=- 1 
ENDIF 

101  CONTINUE 
ENDIF 

50  RETURN 
END 


|-Wr*';V'i-r-,f**-,'f*Vr-.V***********-,V********-'f***Vr*:?fVrVr-.V*********,!V*'!,.***Vr****Vf*'!'r***** 


SUBROUTINE  OUTPUT( S , GAMMA 1 , GAMMA 2 , GAMMA3 ,R,V,B,C,K,CW,D, 
+  BAUDNO , DHAT , BDTYPE , MAX ) 

FOR  BYTE  ERRORS,  WRITE  PERTINENT  DATA  TO  RSOUTPUT  FILE... 


INTEGER  GAMMA 1 , GAMMA2 , GAMMA 3 , B , C , K , MAX 
INTEGER  S(  -  2:  2) ,R( 0:  MAX) ,V(0:  MAX) ,CW(0:  MAX) 

INTEGER  DC 0:  MAX) , BAUDNO, DHAT( 0: MAX) , BDTYPE 

WRITEC  6 ,*) ' . 

+  - - - * 

WRITE (6,*) ' BAUDTYPE= ' , BDTYPE 
WRITE( 6 , * ) ' CODEWORD  NUMBER= ' , BAUDNO 
WRITE( 6 ,904) 

904  FORMAT ( ' SYNDROMES ' ) 

WRITEC 6,905)  S( -2) ,S( -1) ,S(0)  ,S( 1)  ,S(2) 

905  F0RMAT( IX, 12 . 3X, I2,3X, 12 , 3X, 12 , 3X, 12) 

WRITEC 6,*)' 

WRITEC 6,*)'  ' 

WRITEC  6. 906) 

906  FORMATC  GAMMAS') 

WRITEC  6 , 907 )  GAMMA 1 , GAMMA2 , GAMMA 3 

907  FORMATC 5X , 12 , 5X , 12 , 5X, 12) 

WRITEC  6,*) ' B ,C ,K' 

WRITE(6,*)B,C,K 

WRITEC 6,*)'  ' 

WRITEC 6, 908) 

908  FORMATC 1 IX, ' I' ,10X,'DATA' ,7X, ’ENCODE' ,7X,'RECD' ,8X,'ESTCW' , 

+  5X, ' ESTDATA' ) 

DO  71  J=0,MAX-1 

WRITE  C 6 , * )  J,D(J),CW(J),RCJ),V(J), 

+  '  DHAT(J) 

71  CONTINUE 

RETURN- 
END 


Q********V:********************Vf**Vr*Vf**';V***Vf*********Vf****Vf**ilf!V'!V**'!'.:,"'iV** 
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SUBROUTINE  INTRAN( D 1 ,  D ,  MAX ,  NLESSK ) 


SUBROUTINE  TO  TRANSLATE  REAL  DATA  INTO  INTEGERS.  .  . 
(SINCE  CODER/DECODER  DEALS  ONLY  WITH  INTEGERS) 


DO . INTEGER  OUTPUT 

D1 . REAL  INPUT 


INTEGER  MAX, NLESSK 
INTEGER  D(0:  MAX-NLESSK-1) 

REAL  D1(0: MAX-NLESSK-1) 

DO  1  1=0, MAX-NLESSK-1 
DO  2  J=0 ,MAX- 1 

IF  (Dl(I).EQ.  (J/l.  0))  THEN 
D(I)=J 
END  IF 

2  CONTINUE 

1  CONTINUE 
RETURN 
END 


SUBROUTINE  B ITREV(N , M , XTMP , X) 


C  THIS  SUBROUTINE  ACCEPTS  AS  INPUT  THE  COMPLEX  ARRAY 
C  CONTAINING  THE  VALUES  X'TMPQ  WHICH  WERE  READ  FROM 
C  THE  INPUT  FILE.  THE  OUTPUT  OF  THIS  SUBROUTINE  IS 
C  THE  COMPLEX  ARRAY  X( )  WHICH  CONTAINS  THE  INPUT 
C  VALUES  IN  ’BIT-REVERSED’  ORDER. 

C  (FROM  THE  LIBRARY  OF  PROFESSOR  STRUM  AT  NPS) 

C 

C . 


COMPLEX  XTMP(0:  N-l) ,X(0:  N-l) 

DO  10  K=0 ,N-1 
NEWADR=0 
MADDR=K 
DO  20  1=0, M-l 

LRMNDR=MOD( MADDR , 2 ) 

NEW ADR=NE W ADR+LRMN  D  R* 2** ( M - 1 - 1 ) 
MADDR=MADDR/2 
20  CONTINUE 

X( NEWADR )=XTMP( K ) 

10  CONTINUE 


RETURN 

END 
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SUBROUTINE  FFT(N,M,X) 
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THIS  SUBROUTINE  ACCEPTS  AS  INPUT  THE  COMPLEX  ARRAY  X( ) , 
CALCULATES  THE  FAST  FOURIER  TRANSFORM  ( FFT)  OF  THE 
ARRAY,  AND  RETURNS  THE  RESULTING  SEQUENCE  IN  THE 
ORIGINAL  ARRAY  X(). 

(FROM  THE  LIBRARY  OF  PROFESSOR  STRUM  AT  NPS) 


COMPLEX  X( 0 :  N-l) ,W,TMP 
PI =4.  0*ATAN(  1.  0) 

EN=N 


DO  50  L=1 ,M 
ISPACE=2**L 
S=N/ ISPACE 
IWIDTH=I SPACE/2 
DO  40  J=0 , ( IWIDTH- 1) 

R_Q*  J 

ALPHA=2.  0*PI*R/EN 
W=CMPLX( C0S( ALPHA) , -SIN( ALPHA) ) 
DO  30  IT0P=J,N-2, ISPACE 
IBOT=ITOP+IWIDTH 
TMP=Xf IBOT)*W 
X(  IBOT)=Xf ITOP) -TMP 
X(ITOP)=X(ITOP)+TMP 
0  CONTINUE 

0  CONTINUE 

0  CONTINUE 

RETURN 

END 


SUBROUTINE  INVFFT( N , M , X) 

THIS  SUBROUTINE  ACCEPTS  AS  INPUT  THE  COMPLEX  ARRAY 
X()  AND  RETURNS  THE  INVERSE  FFT  OF  THE  ARRAY. 

(FROM  THE  LIBRARY  OF  PROFESSOR  STRUM  AT  NPS) 


COMPLEX  X(0: N-l) 

EN=N 

C  CALCULATE  THE  COMPLEX  CONJUGATE  OF  THE  INPUT  DATA. 

DO  70  1=0, N-l 
X( I )=CONJG(X( I ) ) 

70  CONTINUE 

C  CALCULATE  THE  FAST  FOURIER  TRANSFORM  OF  THE  ARRAY. 

CALL  FFT(N,M,X) 

C  CALCULATE  THE  COMPLEX  CONJUGATE  OF  THE  RESULTING  ARRAY. 
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DO  80  1=0, N-l 

X( I )=CONJG(X( I ) )/EN 
80  CONTINUE 


RETURN 

END 


■•}Wr‘>Wr*Vf'}W<,*Vr  .f>Wryc,>Wr5V,;WnWr*,sV,3V*V?*,}VVr7V,5V,******,jW?*******'jWr,>V*,3VVr3V,>Wc* 


SUBROUTINE  GAUSS( AVG , VAR , PTS , AWGN.TEMP 1 , TEMP2 ) 

GAUSS -SUBROUTINE  TO  OUTPUT  A  GAUSSIAN  RANDOM  VARIABLE  ARRAY  WITH 
MEAN=AVG  AND  VARIANCE=VAR 

AVG . MEAN  OF  THE  GAUSSIAN  RANDOM  VARIABLE 

VAR . VARIANCE  OF  THE  GAUSSIAN  RANDOM  VARIABLE 


INTEGER  PTS 

REAL  AVG , VAR , TEMP 1(0:  PTS-1) ,TEMP2(0: PTS-1) 

REAL  A,B ,PI 

REAL  AWGN(0:  PTS-1) 


PI =4. 0"ATAN( 1.  0) 

CALL  RNUN( PTS, TEMPI) 

CALL  RNUN ( PTS , TEMP2 ) 

DO  10  1=0, PTS-1 

A=-2. 0*ALOG(  TEMP 1(1)) 

B=COS( 2*PI*TEMP2( I ) ) 

AWGN( I ) =SQRT( A ) *B*SQRT( VAR) +AVG 
10  CONTINUE 
RETURN 
END 
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